Non-congruence of the nuclear liquid-gas and deconfinement phase transitions 
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We investigate non-congruent first-order phase transitions (PTs) relevant for heavy-ion collisions 
and neutron stars in Coulomb-less models. Two different phase transitions are considered: the 
nuclear liquid-gas PT at sub-saturation densities and the deconfinement FT at high densities and/or 
temperatures. For the first FT, we use the FSUgold relativistic mean-field model and for the second 
one the relativistic chiral SU(3) model. The chiral SU(3) model is one of the few models for the 
deconfinement phase transition, which contains quarks and hadrons in arbitrary proportions (i.e. 
a "solution") and gives a continuous transition from pure hadronic to pure quark matter above a 
critical point. The study shows the universality of the applied concept of non- congruence for the two 
phase transitions with an upper critical point, and illustrates the different typical scales involved. 
In addition, we find a principle difference between the liquid-gas and the deconfinement FTs: in 
contrast to the ordinary Van-der-Waals-like FT, the phase coexistence line of the deconfinement FT 
has a negative slope in the pressure-temperature plane. Furthermore, we also find that the non- 
congruent features of the deconfinement FT become vanishingly small around the critical point. 

PACS numbers: 05.70.Fh 25.75.Nq 21.65.-f 26.60.-c 
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I. INTRODUCTION 

Nuclear matter is expected to undergo two different 
major phase transitions (PTs) : the liquid-gas phase tran- 
sition (LGPT) of nuclear matter at sub-saturation den- 
sities and moderate temperatures and the deconfinement 
and chiral symmetry restoration phase transition at high 
densities and/or temperatures. For convenience we will 
call the latter also the quark hadron phase transition 
(QHPT) or QCD phase transition. These two PTs are 
actively discussed in the contexts of heavy-ion collisions 
and astrophysics. The latter includes the interior of com- 
pact stars, i.e. neutron stars (NS) or so-called hybrid 
stars which have quark matter in their core. 

Various effective models for nuclear matter which are 
constrained by properties of nuclei have shown that the 
LGPT of bulk uniform nucleonic matter, i.e. consisting 
of neutrons and protons without Coulomb interactions, 
is of first order. A prominent example is Chiral effective 
field theory, which only very recently has been applied to 
investigate the LGPT in asymmetric nuclear matter at 
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sub-saturation densities [l[. Nucleonic matter can also 
be calculated with realistic nucleon-nucleon interactions 
using Monte-Carlo algorithms and sophisticated many- 
body techniques, which have been used to determine the 
LGPT as well, see e.g. Ref. From such theoretical 
studies one can conclude that there is enough reason 
to expect the first-order nature of the LGPT. Further- 
more, there is also experimental evidence for this from 
intermediate-energy heavy- ion collisions . 

A direct solution of the QCD equations for the high- 
density low-temperature region of the phase diagram 
around the deconfinement and chiral symmetry restora- 
tion phase transition is currently out of reach. Ab-initio 
solutions of QCD are not possible, because one is in the 
non-perturbative regime, respectively, they exist only for 
very high densities and/or temperatures j6l-[l0j. Addi- 
tionally, simulations on the lattice have inherent prob- 
lems at finite density. As a consequence, effective models 
for QCD matter have to be used, resulting in different 
varieties of possible QCD phase diagrams [ll| - [27l |. Many 
of these models predict that the QCD phase transition 
at low temperatures is of first order like the LGPT, but 
some also predict a cross-over transition in this regime. 

Phase diagrams of isospin asymmetric matter are of 
extreme importance for the complete understanding of 
QCD and nuclear matter. They arc highly related to 
the symmetry energy, as explained, e.g., in Refs. (28l - 
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1301 1 ■ Such studies are also used to analyze the effect of 
model parameters on the QCD phase diagram [3ll - l33| . 
The effect of different isospin/charge assumptions has 
been studied already extensively in the literature for the 
LGPT [33, Isil - ISTl . and also experimentally [H, [s^ , and 
for the QHPT [3ll^, SO] as well. 

Some authors (H, [H liol-li^ have stated that the 
LGPT and QHPT changes from first to second order (ac- 
cording to the Ehrenfest classification) if one goes to an 
asymmetric system. This was concluded from the non- 
standard behavior of thermodynamic quantities during 
an isothermal crossing of the two-phase region. One of 
the main statements of the present paper is that this 
non-standard behavior in asymmetric matter has instead 
its natural explanation as the typical manifestation of a 
so-called "non-congruent" first-order PT. 

Based on the assumption that both the LGPT and 
the QHPT are of first order, in the present investigation 
we concentrate on the detailed thermodynamic aspects 
of the two phase transitions and investigate their non- 
congruent features. The two PTs are studied for the sce- 
narios of heavy-ion collisions of symmetric and asymmet- 
ric nuclei. For the QHPT, we are also investigating the 
scenario of the interior of a neutron, respectively hybrid 
star. For this purpose we use a relativistic chiral SU(3) ef- 
fective model. This model predicts that both the LG and 
the QCD PT at low temperatures are of first order. Due 
to technical reasons, the chiral SU(3) model is applied 
only for the QHPT. For the description of the LGPT we 
apply the FSUgold relativistic mean-field model. Even 
using only one selected theoretical model for the QHPT 
and one for the LGPT, our main conclusions arc to some 
extent model- independent, because the applied thermo- 
dynamic concepts are rather universal. 

The structure of this article is as follows: in Sec. HIl we 
discuss (non-) congruence of PTs. In Sec. IHII we describe 
the effective model used for the calculations of the LGPT, 
the FSUgold equation of state (EOS). In Sec.|lV]wc con- 
tinue with the description of the chiral SU(3) model for 
the QHPT. In Sec. |V] we specify our thermodynamic 
model and setup used for the two phase transitions and 
the different physical systems. In Sec. I VII we analyze and 
compare in detail the results for the different scenarios, 
with a focus on the structure of the resulting phase di- 
agrams and the non-congruent features. In Sec. IVIII we 
summarize our main findings and draw conclusions. 



II. CONGRUENCE/NON-CONGRUENCE OF 
PHASE TRANSITIONS 



A. Coulomb interactions 

Before discussing non-congruence, it should be stressed 
that for all thermodynamic systems in the present paper, 
including those corresponding to matter in neutron stars. 
Coulomb interactions are not taken into account explic- 
itly, in spite of the presence of charged species (protons. 



quarks, leptons, etc.). This is what we call a "Coulomb- 
less" model description. Also, surface effects are ne- 
glected in our work. As a consequence, the two-phase 
mixtures at equilibrium (i.e., not metastable) within the 
two-phase regions are always described as coexistence of 
two macroscopic phases. 

The simplification of Coulomb-less is to some extent 
reasonable for the theoretical description of relativistic 
heavy-ion collisions, where one has a net electric charge 
but Coulomb energies are small compared to the typi- 
cal collision energies. Furthermore, the long-range na- 
ture of Coulomb forces could be ignored in view of the 
small size of the ensemble of heavy-ion collisions prod- 
ucts. However, for the same reason it is questionable 
whether the thermodynamic limit is fulfilled or not [43j . 
On the other hand, for the description of nuclears clus- 
ters appearing in the nuclear liquid-gas PT of low-energy 
heavy-ion coUsions, Coulomb and surface energies are in 
fact crucial. Nevertheless, the bulk Coulomb-less treat- 
ment gives useful insight into the main characteristics of 
the phase transition. 

Matter in neutron stars has to be overall charge 
neutral in order to be gravitationally bound. In this 
case. Coulomb interactions and corresponding surface 
effects can be included in a more detailed mesoscopic 
description, leading to structured mixed phases. Usu- 
ally these are called the "pasta phases" or the "pasta 
plasma" (e.g. Ref. [iJl) where in the simplest case one as- 
sumes a certain fixed dimensionality/geometry of the two 
phases. In a more sophisticated description, with charge 
screening effects included, one obtains more gen eral non- 
uniform particle distributions (e.g. Refs. j45l - l48l |). Also a 
three-dimensional Hartree-Fock treatment leads to pasta 
phases for the LGPT close to normal nuclear matter den- 
sity 0. 

Without such a detailed treatment, the effect of 
Coulomb interactions in NSs can be estimated by differ- 
ent assumptions for charge neutrality (sol - fs^ . which we 
will use in the present study. In this case, there are two 
limiting cases of the mesoscopic description, depending 
on the model-dependent surface tension (e.g. Ref. [27| ) 
between the two phases, which can be ap pro ximated 
within a macroscopic description of the PT |52| - f54| . In 
the first case one assumes an infinitely high surface ten- 
sion. Thus the two phases have to become infinitely 
large and can be treated as being macroscopic. The in- 
clusion of Coulomb forces then automatically leads to 
local charge neutrality (LCN). In the astrophysics com- 
munity, the PT in beta-equilibrated NS matter of macro- 
scopic phases with assumed local charge neutrality is typ- 
ically called a "Maxwell-PT" . In terrestrial plasmas this 
is denoted more accurately by the Gibbs-Guggenheim 
conditions for phase equilibrium, see e.g. Ref. [55|. In 
the second limiting case of the mesoscopic description 
one assumes a vanishing surface tension. This leads to 
a highly dispersed charged and non-soluble mixture of 
micro-fragments of one phase into the other, a mixed 
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phase, which also could be called a charged "emulsion" 
Due to the vanishing size of the fragments, the Coulomb 
interactions are negligible in this case. If the remaining 
Coulomb and surface effects are neglected, this situation 
can be approximated by two macroscopic phases in phase 
equilibrium which have opposite electric charge and the 
system is only globally charge neutral (GCN). However, 
the thermodynamic and Coulomb-less treatment of such 
micro-fragments remains an approximation. In astro- 
physics, the Coulomb-less PT in beta-equilibrated NS 
matter with GCN and macroscopic phases is typically 
called a "Gibbs-PT" [5^. 

We remark that very often in the astrophysical litera- 
ture, matter in the two-phase coexistence region of any 
PT, including those in neutral systems, is generally said 
to be in a "mixed phase" . We think it is more accurate 
to denote this as a "two-phase mixture" and to reserve 
the term "mixed phase" only for the state of matter ob- 
tained in the mesoscopic description of PTs in Coulomb 
systems with a low surface tension, as described above. 

B. Definition of non-congruence 

The term "non-congruent" (or "incongruent" ) phase 
transition (NCPT) denotes the situation of phase coex- 
istence of two (or more) macroscopic phases with dif- 
ferent chemical compositions (see the lUPAC definition 
[stI ] and Ref. [HH). Such systems are also called "bi- 
nary", "ternary", etc., in contrast to "unary" systems. 
NCPTs are well known since long ago in many terrestrial 
applications as a particular type of PTs (regardless of 
the term), e.g. in low-temperature solution theory (see 
e.g. Ref. [58|), in the theory of simple liquid mixtures 
of hydrocarbons (see e.g. Ref. [11]), or in the theory of 
crystal-fluid and crystal-crystal phase diagrams in chemi- 
cal compounds. NCPTs are also known in nuclear physics 
[sif , in heavy-i on p hysics jooj , and also in the physics of 
compact stars [56[ since quite some time, but the term 
"non-congruent" has been introduced to these areas of 
physics only recently (see below). 

The variants of terrestrial NCPT which arc the closest 
ones to LGPT and QHPT discussed in the present arti- 
cle, are PTs in high-temperature, chemically reacting and 
partially ionized plasmas — typical products of extremely 
heated chemical compounds. NCPTs were studied thor- 
oughly for high-temperature uranium-oxygen systems 
during hypothetical "severe accidents" in the framework 
of nuclear reactor safety problems [6ll - l66| . The univer- 
sal nature of this type of PT and its applicability for 
most astrophysical objects was claimed and illustrated 
in Ref. [6^, using the examples of (hypothetical) plasma 
phase transitions in the interiors of Jupiter and Saturn, 



Another term (culinary like "pasta", "spaghetti", etc.) was pro- 
posed for this emulsion-like mixture: "milk phase" i.e. highly 
dispersed mixture of oil micro-drops in water [55| . 



brown dwarfs, and extrasolar planets. The identification 
that most PTs in neutron stars are non-congruent, in par- 
ticular for the QHPT in hybrid stars, was claimed first 
at several conferences by I.I. and then published recently 
in Ref. Nowadays, the term "non-congruent" phase 
transition is already used in the astrophysical literature 
[41I Our theoretical description of the LGPT and 
QHPT as non-congruent phase transitions in the present 
study is based essentially on experience from terrestrial 
applications. 

Regarding the above standard terrestrial definition of 
NCPTs, it should be noted that coexisting phases could 
have differences in ionization composition besides the dif- 
ference in chemical composition, i.e., atoms or molecules 
could appear in different ionized states. However, this ad- 
ditional degree of freedom does not count in the definition 
of a non-congruent PT when one deals with phase coex- 
istence of two electroneutral macroscopic phases (or a 
mixture of electroneutral macroscopic fragments of both 
phases) and if the Coulomb interactions between charged 
species are included. This exception is due to the va- 
Hdity of LCN in this described above. This is 

different in the framework of the Coulomb-less approach, 
where the PT may be non-congruent, even for a plasma 
of one chemical element. Hence, the definition of "non- 
congruence" in the case of chemically reacting macro- 
scopic Coulomb systems does not depend on the ioniza- 
tion. This is different in the Coulomb-less approach, in 
which positive and negative charges (say nuclei and elec- 
trons) play the role of different chemical elements. 

The generalization of this definition for first-order 
phase transitions in dense nuclear matter, described as 
Coulomb-less systems with more than one conserved 
"charge" , is obvious: phase coexistence of two (or more) 
macroscopic phases with different composition of the 
charges, including electric charge. In the Coulomb-less 
description of nuclear matter, the electric charge is equiv- 
alent to an additional chemical element in the case of 
terrestrial chemical thermodynamics. Conversely, if the 
long-range Coulomb forces are taken into account in a 
macroscopic description (leading to local charge neutral- 
ity), the electric charge is excluded from the number of 
equivalents of chemical elements in the generalized defi- 
nition of NCPTs in dense nuclear matter. In conclusion, 
it is the number of globally conserved charges which is 
relevant for the non-congruence. 

The classification with respect to non-congruence in 
the case of the widely accepted hypothetical mesoscopic 
scenarios in Coulomb systems of nuclear matter (simple 
mixed phase, more complicated structured mixed phase, 
respectively pasta plasma or pasta phases) is still an open 
question. In a strict thermodynamic sense, the state 
of matter in such mesoscopic calculations should not be 
seen as the two-phase coexistence of a first-order PT, but 
rather as a single phase with a non-uniform substructure. 

There is a famous example from the context of neutron 
stars which illustrates the definition of non-congruence: 
in beta-equilibrated, cold neutron stars baryon number 



4 



and total net electric charge (which has to be zero) are 
two conserved charges. As described in detail in Sec. Ill A| 
there are two typical choices for the treatment of charge 
neutrality for PTs of macroscopic phases. In the first 
case, one assumes LCN, and obtains a congruent phase 
transition of a unary system. In astrophysics this is usu- 
ally called the "Maxwell-PT" , which is then used as a 
synonym for congruent phase transitions in general. In 
the second case one assumes GCN and the Coulombless 
approximation. In this case the system is binary and the 
PT is non-congruent [s^. In astrophysics this is often 
called the "Gibbs-PT", and again taken as a synonym 
for non-congruent phase transitions in general. The clas- 
sification with respect to "Gibbs" or "Maxwell" of mat- 
ter in supernovae or proto-neutron stars with possibly 
trapped neutrinos was given in Ref. . Nuclear matter 
in heavy-ion collisions also has more than one conserved 
charge, namely net baryon number, net electric charge 
and also net weak flavor, respectively isospin, due to 
the fast timescales involved. Thus Coulomb-less PTs in 
heavy- ion collisions will in general also be non-congruent, 
see also Ref. [G^I ■ Ref. addresses experimental conse- 
quences of the QHPT as a non-congruent PT. The previ- 
ous arguments are valid for both phase transitions, LGPT 
and QHPT, just the typical scales involved and the quan- 
titative behavior is different. 



C. Characteristics of non-congruent PT 

It was shown in Refs. [b^. IssI. [sgI. l60l - [66| and many oth- 
ers, and it will also be shown below, that non-congruency 
significantly changes the properties of all phase tran- 
sitions, namely: (A) Significant impact on the phase 
transformation dynamics, i.e., a strong dependence of 
the phase transition parameters on the rapidity of the 
transition [63|. (B) The thermodynamics of phase tran- 
sitions becomes more complicated. The essential changes 
include: (1) significant change in properties of the singu- 
lar points (critical point first of all) and separation of 
critical point and endpoints, such as temperature end- 
point, pressure cndpoint, etc. (2) significant change in 
the scale of two-phase boundaries in extensive thermody- 
namic variables (say P-V, V-T, H-T, etc.) and even in 
topology of all two-phase boundaries in the space of all 
intensive thermodynamic variables, i.e., pressure, tem- 
perature, specific Gibbs free energy etc. Note, that this 
is valid for both types of PTs: with and without a crit- 
ical point (e.g. gas-liquid-like PT and crystal-fiuid-like 
PT, correspondingly). One of the most remarkable con- 
sequences of the non-congruence in NCPT is the appear- 
ance of a two-dimensional "banana-like" region instead of 
the well-known one-dimensional P-T saturation curve for 
ordinary (congruent) PTs (see Fig. 1 in Ref. [III). The 
same should be expected in the plane of the widely used 
pair of variables: temperature - baryon chemical poten- 
tial (see below). (3) Closely connected to this is the sig- 
nificant change of the behavior in the two-phase region: 



i.e. isothermal and isobaric crossings of the two-phase 
region do not longer coincide. The isothermal NCPT 
starts and finishes at different pressures, while the iso- 
baric NCPT starts and finishes at different temperatures. 
Basically, the pressure on an isotherm is monotonically 
increasing with density. 

Aspect (3) of non-congruent phase transitions is well 
studied in the context of neutron stars [s^ . Inside a neu- 
tron star, the pressure has to decrease monotonically with 
the radius. A congruent phase transition leads there- 
fore to a spatial separation of the two coexisting phases, 
with a discontinuous jump in density and all extensive 
thermodynamic variables at the transition radius inside 
the neutron star. Conversely, for a non-congruent phase 
transition a spatially extended two-phase coexistence re- 
gion is present, with a continuous change of total density, 
total energy density, etc, throughout. We remark that for 
the LGPT there exist several works which also have dis- 
cussed the other characteristic features of non-congruent 
phase transitions and only used a partially different ter- 
minology, see, e.g., Refs. [33,[3^[37|. For the QHPT, the 
possible non-congruence has not been discussed in such 
detail as is done here. 

Furthermore, even nowadays publications are still ap- 
pearing which do not treat the thermodynamic aspects 
of non-unary phase transitions, i.e., the non-congruent 
features, in a proper way. For example aspect (3) some- 
times lead to the conclusion that one has a second-order 
phase transition according to the Ehrenfest classification, 
as e.g., in Refs. [ssl. [stI. |40| - |42| . However, the two phases 
have different order parameters like densities, entropies, 
asymmetries, etc., and, most importantly for our pur- 
poses, different generalized "chemical" compositions. At 
the interface between the two macroscopic phases there 
is a discontinuous jump of the order parameter and thus 
the phase transition is still of first order. Also the first 
two aspects of (B) from above are sometimes overseen 
or neglected in the literature which means that the non- 
congruence is not fully taken into account (compare, e.g., 
with Ref. [H). 



D. Isospin symmetry, azeotrope 

The isospin symmetry of strong forces plays an impor- 
tant role for the possible non-congruence of the LGPT 
and the QHPT. Independently of density and temper- 
ature, isospin symmetric nuclear matter always repre- 
sents the state with the lowest thermodynamic poten- 
tial (neglecting Coulomb interactions and assuming equal 
masses of protons and neutrons). Thus the isospin chem- 
ical potential is zero for symmetric nuclear matter. As a 
consequence, isospin does not appear as a relevant charge 
for symmetric nuclear matter because this degree of free- 
dom is not explored, i.e. even in a first-order phase tran- 
sition the involved phases remain symmetric. Therefore 
the LG and QH phase transitions remain congruent if the 
system is exactly symmetric and if no other globally con- 
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served charge than baryon number is involved (see also 
Appendix [B|. This is called "azeotropic" behavior, de- 
noted for a system with more than one conserved charge 
whose charge ratios cannot be changed by distillation for 
a certain azeotropic composition. The ensemble of such 
azeotropic points in the parameter space, e.g., for all tem- 
peratures, is called an azeotropic curve. Note that the 
isospin asymmetry of the hot state of matter in a heavy- 
ion collision experiment is mainly set by the initial charge 
to mass ratio ZjA of the colliding nuclei, e.g., Zj A ~ 0.4 
in Au-I-Au collisions. 

E. Unified EOS 

As mentioned in point (1) in Section lll C[ another con- 
sequence of non-congruent phase transitions is the possi- 
ble emergence of critical points, which are different from 
the points of maximum temperature, pressure, or ex- 
tremal chemical potential. To obtain such critical points 
and endpoints at all, it is necessary that both of the two 
involved phases are calculated with the same theoretical 
model ("unified" or "single" EOS approach). In other 
words, one has to use only one underlying many-body 
Hamilton operator. This is in contrast to a "two-EOS" 
approach, where two different EOS models are applied 
for the two phases in coexistence. Such a "two-EOS" 
description can have several short-comings, as it cannot 
contain critical points and endpoints (see Appendix [A} 
and standard literature, e.g., Ref. ^M) and it does not 
give a consistent description of meta-stable or unstable 
matter in the binodal, respectively spinodal regions, e.g., 
for a liquid-gas type PT. In summary, in the "unified" 
EOS approach both coexisting phases are presumed as 
isostructural (like gas and liquid) with a possible contin- 
uous transition from one phase into another, while in the 
two-EOS approach this is impossible. 

Almost all studies of the LGPT are based on the uni- 
fied EOS approach. This also applies for our investigation 
of the LGPT with the FSUgold relativistic mean-field 
model. Unified EOS approaches for the QHPT are usu- 
ally either built with only hadrons or only quarks. Thus 
they do not give the expected degrees of freedoms for 
one of the two phases. Alternatively, often the two-EOS 
approach is applied for the QHPT (see, e.g., Ref. fioj ) 
to have the right degrees of freedom. On the other hand 
this approach cannot contain all possible non-congruent 
features of the singular points, as explained above. The 
chiral SU(3) model used in the present study is one of the 
few unified-EOS approaches for the QHPT that contains 
hadronic as well as quark degrees of freedom. These can 
appear, in principle, in arbitrary proportions (solution- 
like mixturqj) with the interactions leading to the correct 



^ Another term was proposed for this solution-hke mixture: 
"vodka phase", i.e., a solution of spirit in water with arbitrary 
proportion [55| . 



behavior for low, respectively high, densities and temper- 
atures. See Refs. [la, [t^I for another unified-EOS model 
that also contains hadronic and quark degrees of freedom. 
As another exception of a unified-EOS approach for the 
QHPT with the correct degrees of freedom there is the 
EOS of Ref. [7l|, where the two-EOS approach is trans- 
formed into a one-EOS version with the use of a special 
spline-based interpolation procedure. 



III. FSUGOLD RMF MODEL 

For the liquid-gas phase transition (LGPT) of nucle- 
onic matter we apply a relativistic mean-field (RMF) 
model. In principle, also the Chiral model could be used 
for this, as it also contains the liquid-gas phase transition 
[t^ . However, due to the different characteristic scales 
involved and for numerical reasons, we use a dedicated 
model for the liquid-gas phase transition which occurs at 
sub-saturation densities. We choose the FSUgold RMF 
parameterization [t^, because of its excellent descrip- 
tion of matter around and below saturation density and 
because its neutron matter EOS is in agreement with 
recent experimental and observational constraints (see 
e.g. Ref. (ZJI)- Its Lagrangian is based on the exchange 
of the isoscalar-scalar ct, the isoscalar-vector uj and the 
isovector- vector p mesons between nucleons. Particular 
for FSUgold, also the coupling between the uj and the p 
meson is included. This leads to a better description of 
nuclear collective modes, the EOS of asymmetric nuclear 
matter, and a different density dependence of the symme- 
try energy [TSj . The free parameters of the Lagrangian, 
the meson masses and their coupling constants, are deter- 
mined by fits to experimental data, more specifically to 
binding energies and charge radii of a selection of magic 
nuclei. 

The only baryonic degrees of freedom in FSUgold are 
neutrons and protons. For the typical densities and tem- 
peratures of the LGPT, hyperonic or quark degrees of 
freedom are not relevant. Because FSUgold considers 
only the "elementary" particles of the liquid-gas phase 
transition but not any compound objects, respectively 
bound complexes like light or heavy nuclei, it belongs to 
the class of so-called physical descriptions of phase tran- 
sitions. In this description, all effects of bound complexes 
are presumed to be taken into account by the interactions 
("non-ideality") of the "elementary" particles (see, e.g., 
Ref. [zl]). However, for a more detailed description of the 
nuclear EOS like, e.g., used in simulations of core-collapse 
supernovae, the formation of nuclei and nuclear clusters 
has to be incorporated explicitly, see e.g., Refs. fi?!. [ttI- 
[83| . For high temperatures, light nuclei like the deuteron 
or alpha particle are most important, whereas at low tem- 
peratures, heavy and also super-heavy nuclei give the 
dominant contribution. If all possible compound ob- 
jects (i.e. nuclei) were included as a chemical mixture, 
one would obtain a quasi-chemical representation, as e.g. 
done in Ref. [s^l • This can cause substantial changes and 
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even to a quenching of the hquid-gas phase transition as 
a first-order phase transition, see also Ref. [s^l- How- 
ever, even in this case one can use the analogy between 
the characteristic changes of the nuclear composition and 
the behavior of the gas and the liquid phases in a pure 
thermodynamic treatment, see e.g. Refs. [t^H^. It is 
confirmed in many studies, that the mean-field without 
clusterization overestimates the region of instability, see, 
e.g., Refs. [7^, 113. Because Coulomb-interactions and 
clusterization are more important in the cold catalyzed 
matter of neutron stars than in the hot plasma of heavy- 
ion collisions, we are discussing the LGPT only in the 
latter scenario. 



IV. CHIRAL SU(3) MODEL 

The non-linear realization of the sigma model [s^, 
is built on the original linear sigma model [89l.l90|. includ- 
ing the pseudo-scalar mesons as the angular parameters 
for the chiral transformation, to be in better agreement 
with nuclear physics results. It is an effective quantum 
relativistic model that describes hadrons interacting via 
meson exchange, similar to the FSUgold RMF interac- 
tions. However, the model is constructed in a chirally 
invariant manner as the particle masses originate from 
interactions with the medium and, therefore, go to zero 
at high densities/temperatures. 

The Lagrangian density of the model in the mean- 
field approximation (all particles contribute to the global 
mean-field interactions and are in turn affected by them), 
constrained further by astrophysical data, can be found 
in Refs. iQli - iQSj . In this work, we are going to use an ex- 
tension of this model called the Chiral SU(3) model, that 
also includes quarks [72|. Up, down, and strange quarks 
and the whole baryon octet are considered. In this case 
the populated degrees of freedom change from hadrons 
to quarks and vice- versa through the introduction of an 
extra field $ in the effective masses of the baryons and 
quarks. The scalar field $ is named in analogy to the 
Polyakov loop [Q^I since it also works as the order param- 
eter for deconfinement. In this way, the particle masses 
are generated by the scalar mesons whose mean-field val- 
ues correspond to the isoscalar-scalar (a) and isovector- 
scalar {6) light quark-antiquark condensates as well as 
the strange quark-antiquark condensate (C). In addition, 
there is a small explicit mass term Mq and the term con- 
taining 

Mb = gBcrO- + gssTsS + gscC + + cjb-s-^^, (1) 

= gga(T + gqsTsS + g,^C + Mq^ + .9g*(l - (2) 

We remark that for FSUgold only the term with the 
sigma field (with a minus sign) and a large explicit mass 
term Mq^ equal to the nucleon vacuum mass, would 
be present in Eq. ([T|). For FSUgold, the contribution 
of the sigma field is zero in the vacuum and decreases 



the effective mass for finite density. In the Chiral SU(3) 
model, the explicit mass term is much smaller, and the 
nucleon mass in the vacuum is generated mainly by the a 
field (non-strange chiral condensate). With the increase 
of density/temperature, the a field (non-strange chiral 
condensate) decreases from its high value at zero den- 
sity, causing the effective masses of the particles to de- 
crease towards chiral symmetry restoration. The field $ 
assumes non-zero values with the increase of tempera- 
ture/density and, due to its presence in the baryons' ef- 
fective masses, suppresses their appearance. On the other 
hand, the presence of the <& field in the effective mass of 
the quarks, included with a negative sign, ensures that 
they will not be present at low temperatures/densities. 
The potential for $ reads: 

U = (aoT^ + aiti% + azTVI)*' (3) 
-haaTg^ log (1 - + 8$3 _ 3^4^ 

It was modified from its original form in the PNJL model 
[95I [96| in order to also be used to study low temperature 
and high density environments (besides high tempera- 
ture and low density environments). It is a simple form 
to extend the original potential to be able to reproduce 
the physics of the whole phase diagram. Because U now 
also depends on the baryon chemical potential /^s, it will 
provide an extra contribution to the total baryon den- 
sity. It was shown in Refs. [l^,!!^,!!^! that our choice for 
the potential J7($) can also be used in the PNJL model, 
successfully reproducing QCD features. Note that the 
finite-temperature calculations include the heat bath of 
hadronic and quark quasiparticles and their antiparticles 
within the grand canonical potential of the system. 

The coupling constants of Eqs. (I!])-© can be found in 
Ref. [t^I. They were chosen to reproduce lattice data 
as well as information about the phase diagram from 
Refs. [il,!!!,!!!,!!!]. At the critical endpoint a first-order 
phase transition lines begins. The line is calibrated to ter- 
minate on the zero temperature axis at four times satura- 
tion density for charge-neutral beta-equilibrated matter. 
In this way we can reproduce a hybrid star containing 
a quark core. The behavior of the order parameters and 
the resulting phase diagrams will be discussed in Scc. lVIl 

The most important aspect of the chiral SU(3) model 
is that hadrons are included as quasi-particle degrees of 
freedom in a chemical equilibrium mixture with quarks. 
Therefore, the model gives a quasi-chemical represen- 
tation of the deconfinement phase transition (so-called 
"chemical picture" in ter ms o f electromagnetic non-ideal 
plasmas, see, e.g., Ref. jlOOf ). As explained in the in- 
troduction, it is very important for our study that this 
model contains the right degrees of freedom of low and 
high densities, namely hadrons and, respectively, quarks, 
in arbitrary proportions and gives at the same time the 
deconfinement phase transition in a "unified EOS" or 
"single EOS" description. 
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TABLE I. Definitions of total net quantum numbers or to- 
tal net charges which are possibly conserved, depending on 
the scenario considered, and the corresponding chemical po- 
tentials. F = F{T, V, B, S, Q) denotes the total free energy 
which is a function of temperature T, volume V and the total 
net charges B, S, and Q. 



quantity 



definition chem. potential 



baryon number 


B = 


x;, Aj6j 


dF 

as 


T,V,S,Q 


strangeness 


S = 




dF 

^ as 


T,V,B,Q 


electric charge 


Q = 




aF 
^ ag 


T.V.B.S 


electric charge fraction 


Yq 


= Q/B 


not used 



TABLE II. Definitions of net quantum numbers or net charges 
and corresponding chemical potentials of the individual 
phases inside the phase coexistence region. The free energy 
F' of phase I is understood as F' = F' [T' , V' ,B' ,S' ,Q') 
(analogous definitions for phase //). 



quantity 


definition 




chem 


pot. 


baryon number 






dF' 

'aW 


Tyi ,S' ,Q' 


strangeness 


s' = E. K's^ 


I 

Ms 


bf' 
^ as' 


T,V' ,B' ,QI 


electric charge 




Mq 


dF' 

^ aq' 


T,V' ,B' ,S' 



electric charge fraction Yq = /B' 



not used 



V. THERMODYNAMIC SETUP 
A. Definitions 

For the (Coulomb-less) scenarios we are interested in, 
the following three quantum numbers of each particle 
species i are relevant: baryon number 6^, electric charge 
number qi and strangeness Si . The corresponding value s 
can be found in standard textbooks, or e.g. in Ref. |lOl| . 
The quantum numbers of each particle species i also set 
the total net quantum numbers or total net charges of 
the thermodynamic system if the total net numbers of 
particles Ni of each species i are known. The total net 
number Ni is the difference between the number of parti- 
cles and the number of corresponding anti- particles of the 
whole system. The possibly conserved total net quantum 
numbers (which arc extensive) are listed in Table IH Very 
often instead of the electric charge number Q, the inten- 
sive cliarge-to-baryon ratio is used, which is defined in 
the last row of the table. For each of the extensive quan- 
tum numbers a corresponding chemical potential can be 
defined. These are listed in the third column of Table HI 
Later we will also use the following chemical potential fl, 

dF 

= A*i3 + Yq^iQ , (5) 

which is equal to the Gibbs free energy per baryon (see 
Appendix [D| . 

For a state which is inside the two-phase coexistence 
region, two spatially separated macroscopic phases are 
present. Each phase has its own set of extensive ther- 
modynamic variables and chemical potentials, listed in 
Table HIl The total extensive quantities F,V, B, S,Q, Ni 
are given as the linear sums of corresponding quantities 
of the coexisting phases. Particle numbers are connected 
to particle number densities through the volumes of each 
phase: 

pi = Ni/V', 
pi^ = Ni^/V'^. (6) 



B. Scenarios and constraints 

Next we are going to define the different cases of the 
two phase transitions studied in different physical sce- 
narios. An overview of these scenarios is given in Ta- 
ble mil We consider phase transitions in three different 
physical systems: the liquid-gas phase transition of nu- 
clear matter (LG), e.g. in low energy heavy- ion collisions, 
the deconfinement phase transition in high energy heavy- 
ion collisions (HI), and the deconfinement phase transi- 
tion in neutron stars (NS). For the first two scenarios 
LG and HI we investigate symmetric (S) nuclear mat- 
ter with Yq =0.5, and asymmetric (AS) nuclear matter 
with Yq = 0.3. The two different electric charge frac- 
tions correspond to heavy-ion reactions of nuclei with 
different charge to mass ratios Z/A. For ^^^Au, which 
is commonly used in heavy-ion experiments, one has 
Z/A ~ 0.4. However, for peripheral collisions Yq ~ 0.35 
can be reached at certain stages of the evolution as dis- 
cussed in Ref. [1^ . For all of the asymmetric configura- 
tions wc also include a forced-congruent (fc) variant of 
phase equilibrium [ssl IgiI . Igs} , where the composition of 
all conserved charges is forced to be equal in the coex- 
isting phases in frames of Maxwell conditions. In par- 
ticular, the charge fraction is constrained locally. For 
the (Coulomb-less) scenarios of NSs, we investigate the 
effect of local (NSLCN) and global charge neutrality (NS- 
GCN). Next we explain the physical meaning of all of the 
constraints in more detail. 

We remark again that we consider only coexistence 
of macroscopic phases and that we do not consider any 
Coulomb interactions despite the significant participa- 
tion of electrically charged particles, as discussed in the 
introduction. Nevertheless, the electric charge is an im- 
portant quantity for our investigations because it is one 
of the conserved charges which determine the possible 
non-congruence. Furthermore, the electric charge is also 
related to isospin. Let us assume that also the total net 
baryon number B and the total net strangeness S are 
kept constant, just like in all scenarios of LG and HI. 
The quantum numbers of the baryons are directly given 
by the sum of the quantum numbers of their constituent 
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TABLE III. Constraints and particles used in the six difTerent scenarios. For each particle also the corresponding anti-particle 
is included. If a quantity of Table HIl is not listed, its value is not constrained additionally. Note that the combined conservation 
of B, S and Yq is equivalent to the conservation of baryon number, strangeness and isospin. 



case 








constraints 








considered particles 


LGS 


B 


= const. 






= 




Yq = 0.5 




neutrons, protons 


LGAS 


B 


= const. 


S' 




= 




Yq = 0.3 




neutrons, protons 


LGASJc 


B 


= const. 




^s" 


= 






0.3 


neutrons, protons 


HIS 


B 


— const. 


s' 


= s" 


= 




Yq = 0.5 




baryon octet, quarks 


HIAS 


B 


— const. 




= s" 


= 




Yq = 0.3 




baryon octet, quarks 


HIASJc 


B 


— const. 


s' 


^s" 


= 






0.3 


baryon octet, quarks 


NSLCN 


B 


= const. 












: 


baryon octet, quarks, leptons 


NSGCN 


B 


= const. 










Yq = 




baryon octet, quarks, leptons 



quarks. Therefore the total numbers of u-, d- and s- 
quarks (free or bound in baryons) are fixed by the total 
net baryon number B, strangeness S and electric charge 
number Q. If the latter three quantities are kept con- 
stant, the total quark content does not change, i.e. flavor 
is conserved. This means no weak reactions occur and 
also the total isospin of the system is conserved. 

In heavy-ion reactions, the typical timescales are on 
the order of 10~^^ s which is much less than weak re- 
action timescales. Therefore we do not allow for weak 
reactions in the cases LG and HI. This is implemented 
via a fixed value of Yq, conservation of baryon num- 
ber B and conserved total net strangeness S = i). In 
addition to global conservation of the electric charge in 
LGS, LGAS, HIS, and HIAS, we also consider locally con- 
strained charge fractions in the forced-congruent cases of 
LGAS-fc and HIAS_fc. S is set to zero, because initially 
there is no strangeness in the two colliding nuclei. In 
principle, there is still the possibility, that one has net 
strangeness in the two phases with = ~S^^ which is 
known as strangeness distillation [60j . Here we suppress 
this degree of freedom to avoid a ternary phase transi- 
tioi^l and set = S^^ = for simplicity. For HIS, 
HIAS and HIAS_fc this means that the total number of 
strange quarks (free or bound in baryons) is equal to 
the number of anti-strange quarks locally and that there 
is a non-zero strange chemical potential, with two dif- 
ferent values in the two phases. For LGS, LGAS, and 
LGAS_fc strangeness is not relevant at all, because no 
strange particles arc considered, but only neutrons and 
protons. This is appropriate for the typical low energies 
where the nuclear liquid-gas phase transition is relevant. 
We do not consider leptons in the cases of LG and HI, be- 
cause they are not present in the initial configuration and 
their plasma in the later evolution with equal amounts 
of particles and antiparticles would not affect the equi- 
librium conditions between baryons and quarks. 



^ With 'ternary" we mean that one had three globally conserved 
charges with three chemical equilibrium conditions. 



It is important to note that at high temperature, the 
inclusion of light real mesons, like pions and kaons, in 
the system is of extreme importance for thermodynamic 
quantities, since light particles dominate in such regime. 
However, if their interactions with the baryons are neg- 
ligible and if they are not affected by local constraints, 
the inclusion of mesons has no impact on baryonic phase 
diagrams. Here we are concentrating on the baryonic 
component and a more detailed treatment of mesons will 
be postponed to future work. 

In cold neutron stars one typically assumes that all 
possible reactions have reached full equilibrium. Weak 
reactions do not conserve strangeness and, therefore, it 
is not listed as a conserved quantity in Table IIIII for the 
two cases of neutron stars, (NSLCN and NSGCN). In 
principle, weak reactions conserve lepton numbers but in 
cold neutron stars neutrinos can escape freely and, there- 
fore, the interior lepton numbers are also not conserved. 

Finally, electrically charged matter cannot exist in neu- 
tron stars on a macroscopic scale, because otherwise they 
would explode, as Coulomb interactions are many orders 
of magnitude stronger than gravity. Thus we also in- 
clude the lepton contribution in form of electrons and 
muons, which is done easily as they are well described as 
ideal Fermi-Dirac gases. We implement electric charge 
neutrality in two different ways, as discussed in the in- 
troduction. This is done either via enforced local charge 
neutrality (NSLCN), where both macroscopic phases are 
charge neutral and Coulomb forces are absent, or via 
global charge neutrality (NSGCN) in a Coulomb-less de- 
scription, where each of the two phases carries a net elec- 
tric charge which sum up to zero. 

We remark that the scenarios LG and HI described 
above could also be taken as simplified examples of su- 
pernova matter, for which one has similar values of Yq. 
On the other hand, supernova matter has to be charge 
neutral, like matter in neutron stars, and therefore nega- 
tively charged leptons have to be included. For GCN and 
the Coulombless approximation, charged leptons would 
not influence the behavior of the PTs in cases HI and 
LG. However, for a realistic description of the LGPT in 
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supernovae the Coulomblcss bulk treatment is not suf- 
ficient, and the formation of nuclei and nuclear clusters 
has to be taken into account, as noted before. 



C. Phase and chemical equiUbrium conditions 

Based on the previous constraints, the equilibrium con- 
ditions can be derived. First we consider the system out- 
side of the phase coexistence region. If there are more 
particle species than conserved charges, conditions for 
chemical equilibrium are necessary. The chemical poten- 
tial Hi of particle i is related to the chemical potentials 
of the total charges: 



^J■i = hfJ-B + Sifis + qifJ-Q 



(7) 



which allows to calculate the abundances of all parti- 
cles, if the values of the total charges are known. Note 
that for NSLCN and NSGCN, the non-conservation of 
strangeness leads to fis = 0, which is nothing but the 
minimization of the thermodynamic potential with re- 
spect to strangeness. 

We remark that it is also possible to formulate the 
equilibrium conditions of Eq. ([7]) by using chemical po- 
tentials of three selected particles instead of the chemical 
potentials /is, /ig, and fis- We want to give an exam- 
ple for this alternative formulation. Taking the chemical 
potentials of neutrons, protons and lambdas as the basic 
units one obtains from Eq. ([7|: 



(8) 



This sets the chemical potentials of all particles, if /i„. 
Hp, and ha are determined according to the external con- 
straints (see Table HlH) . For example this would lead to: 



He- = 2^A - Hp ■ 



(9) 



For NSs, where leptons are considered, one would also 
get: 



He= Hli= fJ-r. 



(10) 



because of the assumption of non-conservation of the lep- 
ton numbers. 

Inside the phase coexistence region one has to consider 
equilibrium conditions between the two phases. Thermal 
and mechanical equilibrium are given by: 



P = 
T ■ 



pii 



(11) 
(12) 



Inside each phase, one still has relations analogous to 
Eq. (I7|). They give the chemical potential of particle i in 
phase /, respectively //, expressed by the local chemical 
potentials of the charges: 



hI = h,HB 



s,Hs 4 

- SiH^g 



(13) 



Next, one has the chemical equilibrium conditions be- 
tween the two phases. In Coulomb-less sytems, which are 
equivalent to terrestrial chemically reacting systems (e.g., 
Ref. [13), the local chemical potentials of all species in 
coexisting phases must be equal, i.e., hI = I^V i if no local 
constraints are applied, according to the traditional laws 
of chemical thermodynamics. In this case hI = lA^ would 
also follow from /i^ — Hb ■ ~ I^q ■ l^s ~ ' ^^^^ 
Eqs. 113]). However, due to the local constraints applied 
(see Table lIlH) . the inter-phasc chemical equilibrium con- 
ditions depend on the scenario considered, and have to 
be derived, e.g., by means of Lagrange- multipliers (see 
also Ref. [6g). In the following, we list the inter-phase 
chemical equilibrium conditions for the different cases. 

LGS, LGAS, HIS, and HIAS 



Note that Hs order to have 

In the alternative formulation from above, 
(fTS)) would be equivalent to: 



= 5" 

Eqs. (HI 



NSLCN 



Ms = Ms 



Ms 



Ms 







(14) 
(15) 

= 0. 

and 

(16) 
(17) 



(18) 
(19) 



The latter relation comes from the non-conservation of 
strangeness and implies that there is a net strangeness in 
both of the two phases. Note that: 

Mq ^ Mq • 
This means for example: 

Mp 7^ Mp ) 

a4 7^ f^'J ■ 



(20) 

(21) 
(22) 



We remark that according to the Gibbs- Guggenheim con- 
ditions (see for example Ref. [sslp. for a macroscopic 
equilibrium Coulomb syst em o ne should introduce the 
electro- chemical potential jl02l| Hq = I^q + ^^alvani — 
Hq = Hq^ (relative to an arbitrary constant in uniform 
Coulomb systems). With this description, the general- 
ized electro-chemical potentials of all charged particles 
would be equal in the two coexisting macroscopic phases, 
but this is not used here. 
NSGCN 



Ms 



Ms 



Mq = Mq 



So here we have: 



Mp 



h'-p 
- ,,11 



(23) 
(24) 
(25) 

(26) 
(27) 
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LGAS-fc and HIASJc Next, we give the equilibrium 
conditions if the local charge fractions are constrained 
to have the same value, Yq = Yq^ {— Yq). Because 
in the considered cases only baryon number remains as a 
globally conserved charge, the Maxwell construction for a 
congruent phase transition can be used. It is well known, 
that for the "Maxwell" phase transition in a neutron 
star with local charge neutrality and beta-equilibrium the 
baryon chemical potential, which in this case is equivalent 
to the neutron chemical potential, has to be equal in the 
two phases, see Eq. (fT8|) . For HIAS one obtains instead 
the following inter-phase chemical equilibrium condition 
ffii: 



.II 



with the local Gibbs free energy per baryon 



M^ 



dB' 



Ms 



Yr 



Qt^Q 



(28) 
(29) 



(30) 
(31) 



and the analogous expression for fl^^ of phase //. 
Eq. ([^5)1 expresses the equality of the specific Gibbs free 
energy of the two phases, respectively the Gibbs free 
energy per baryon used here (see Appendix |D]). This 
is merely the standard Maxwell construction for a con- 
gruent phase transition, which is also applicable for the 
forced-congruent case. 

In general, the baryon and charge chemical potentials 
will not be the same for the two phases in the phase co- 
existence region, because Eq. (PH)) is the only chemical 
equilibrium condition for cases LGAS_fc and HIAS_fc. 
For a better comparison with the non-congruent variants 
LGAS and HIAS, we will show the phase diagrams of 
LGAS, LGASJc, HIAS, and HIASJc not only as a func- 
tion oi ^B, but also as a function of jl. 

Total chemical potentials inside the two-phase mixture 
The equilibrium conditions given above allow one to de- 
termine the phase boundaries and fully specify the prop- 
erties of the two phases in equilibrium. However, the 
non-equality of local chemical potentials due to local con- 
straints leads to the following complication: in this case 
it is not obvious how the total chemical potentials of the 
charges in the two-phase mixture (defined analogously to 
the ones in Tabled with the local constraints of Table UTTl 
in addition) are related to the local chemical potentials 
of Table [Hi which can have different values in the two 
phases. These relations are derived in Appendix ICl We 
are not aware that these expressions have been published 
in the literature before. 



VI. RESULTS OF CALCULATIONS 

In this section we are showing the results for the phase 
diagrams of each case studied, whereas we begin with the 
LGPT and continue with the QHPT. 




910 912 914 916 918 920 922 924 

/iB [MeV] 

FIG. 1. Phase diagram in the temperature - baryon chemical 
potential plane for case LGS {Yq — 0.5). The saturation 
curve (SC) coincides with the boiling curve (BC). The black 
dot marks the ordinary VdW-like critical point (CP). 



A. Nuclear liquid-gas phase transition 

Fig. [1] shows the phase diagram of case LGS, i.e. 
for the liquid-gas phase transition of symmetric nucle- 
onic matter. In principle, symmetric nuclear matter is 
a two-component, binary system of protons and neu- 
trons, respectively baryon number and isospin. How- 
ever, the nuclear interactions and isospin symmetry lead 
to azeotropic behavior, i.e., the ratio of protons to neu- 
trons does not change during phase coexistence and the 
two coexisting phases remain symmetric. The electric 
charge chemical potential /xg of symmetric nuclear mat- 
ter is zero, independently of density and temperature. 
Therefore no isospin distillation occurs, i.e. there is no 
transfer of isospin per baryon, respectively Yq, between 
the two phases. Since fig = 0, the relation of chemical 
equilibrium with respect to changes of Yq, Eq. psp . is 
automatically fulfilled, and only Eq. (|14p carries relevant 
information. Consequently, symmetric nuclear matter 
behaves like a unary system and the phase transition is 
of congruent type with a phase-coexistence line in the T- 
/is-plane shown in Fig.[TJ This line can be obtained with 
a Maxwell construction by the corresponding constraints 
of Sec. IV C[ Note that the saturation curve (SC) (which is 
also called "dew-point line" ) and the boiling curve (BC) 
(which is also called "bubble-point line" ) coincide in the 
case of congruent phase transitions or azeotropic com- 
positions, and are split into separate boundaries in the 
general case of NCPT. The critical point (CP) marked by 
the black dot, which is also a (critical) endpoint here, is 
located at a temperature of 14.75 MeV and baryon chem- 
ical potential of 912.4 MeV (further values are given in 
Table ITV)) . It is known from other studies, that the CP 
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1 0"^ °1 0"^ 1 0'^ 1 0'^ 1 0"^ 1 0"^ 1 0"^ 1 0'^ 1 0'^ 1 0"^ 1 

P [MeV/fm^] 

FIG. 2. Phase diagram in the temperature - pressure plane 
for case LGS, otherwise the same notation and depiction as 
in Fig. [T] is used. 



of LGS is usually also the global maximum of the phase 
transition temperature, i.e. for all possible values of Yq. 

Ill Fig. [2] we show the pressure-temperature phase dia- 
gram, where we also obtain a phase transition line. Note 
that the pressure on the coexistence line goes to zero in 
the zero-temperature limit. For a congruent phase tran- 
sition the Clapcyron-equation is valid: 

dP - s^^ 

df = W~WS ' ^^^^ 

which describes how the slope of the pressure- 
temperature phase transition line is related to the differ- 
ence in baryon number density and entropy per baryon 
of the two phases. In our investigation, we always have 
Pb ^ Pi/; i-^- first phase is assumed to have lower 
density. In Fig. [2] we see that dP/dT > 0, and thus 
> (where we have replaced "/" by "G" and "//" 
by "L" ) . The gas phase has a higher entropy per baryon 
and is always less dense than the liquid phase, which is 
a characteristic of the liquid-gas phase transition. 

In Fig. [3] the binodal region is shown in the 
temperature-density plane. The gray line to the left of 
the critical point depicts the SC, where droplets of liquid 
form within the nucleon gas. The black line is the BC, 
where bubbles of gas form inside the liquid. The region 
enclosed by the two lines is the phase coexistence region, 
where a two-phase mixture of gas and liquid is present. 
Here, and also in all following plots, filled areas corre- 
spond to states of such a two-phase coexistence. Due to 
the congruent behavior of LGS, for each point inside the 
binodal or phase coexistence region, the gas state on the 
SC is in coexistence with the liquid state on the BC at the 
same common temperature. Thus the gas and the liquid 
arc distinguished from each other by density, whereas the 
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1 0"^ °1 0"^ 1 0'^ 1 0"^ 1 0"^ 1 0"^ 1 0"^ 1 0"^ 1 0"^ 1 0"^ 1 

FIG. 3. The binodal line which encloses the coexistence region 
(filled with gray) in the T - ps plane for case LGS. The black 
dot marks the critical point (CP). The gray line to the left of 
the CP is the saturation curve (SC), and the black line to the 
right is the boiling curve (BC). 



liquid is always more dense. At the critical point the two 
phases are identical. Inside the phase coexistence region, 
the volume fraction of the liquid phase a = /V and 
the gas phase (1 — a) are set by the total baryon number 
density ps = B/V through: 

PS = pi(l - a) +p|a . (33) 

Obviously, one has a = on the SC and a = 1 on the 
BC. 

For asymmetric nuclear matter in the case LGAS one 
obtains a non-congruent phase transition, which can be 
seen in Fig. 31 depicted by the red and blue thick lines. 
For the non-trivial solution of the equilibrium conditions 
in the non-congruent case we have used the method de- 
scribed in Ref. [sHi. The thin gray and black lines show 
the forced congruent variant LGAS_fc which will be dis- 
cussed later. Also in all following plots, we will use col- 
ored solid lines for non-congruent phase transitions. For 
LGAS one has a phase coexistence region in T-fig, en- 
closed by the red and blue thick lines, instead of a single 
line as in Fig. [T]for LGS. Thus one can also distinguish 
between the branch belonging to the SC and the branch 
belonging to the BC. 

As it was stressed in Refs. (ssl. [65|. in non-congruent 
VdW-like phase transitions of gas-liquid type there is no 
more a unique "critical endpoint" . Instead, three sepa- 
rate endpoints exist: maximal temperature (criconden- 
thermj^, maximal pressure (cricondenbar j^, and point 



* The temperature endpoint (TEP) or point of maximal tempera- 
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FIG. 4. Phase diagram in tlie temperature - baryon chemical 
potential plane for the two asymmetric systems LGAS and 
LGASJc with Yq = 0.3. The thin gray and black lines show 
the forced-congruent variant LGAS_fc, and the corresponding 
black dot is the pseudo-critical point. The thick colored lines 
depict case LGAS, where Yq is not constrained locally. The 
corresponding black dot is the critical point and the open 
square the temperature endpoint. The inlay shows a zoom-in 
to the region around the critical point of the non-congruent 
case LGAS. The filled areas are the coexistence regions, (color 
version online) 



of extremal chemical potential. In NCPT, these three 
"topological" endpoints are separated from the singular 
thermodynamic object — the true non-congruent critical 
pointQ. Note also that the critical point of a congruent 
phase transition is determined by: 



dP 



dpi 



dpi 







(34) 



In contrast, for a non-congruent phase transition this cri- 
teria is not applicable, and the critical point does not 
fulfill it in general. 

The inlay of Fig. |4] shows that the temperature end- 
point is different from the critical point. For LGAS, the 



tare, which is also called the "cricondentherm" [59l . 110311 . is de- 
fined as the point with the highest temperature where phase 

coexistence is possible. 

^ The pressure endpoint, also called "cricondenbar" [59l. Il03l |. is 
defined as the point on the binodal where the maximal pressure 
is obtained. 

^ The chemical potential endpoint or point of extremal chemical 
potential is defined as the point where the chemical potential of 
the binodal surface is extremal with respect to temperature. 
The critical point is defined as the point on the binodal surface 
where the two phases are identical. Because it is located on the 
binodal, an infinitesimal change of the state can lead to phase 
separation into two phases which can be distinguished from each 
other by an order parameter. 



FIG. 5. Phase diagram in the plane of temperature and Gibbs 
free energy per baryon jl for case LGAS and LGAS_fc. Oth- 
erwise the same notation and depiction as in Fig. U is used. 



critical point is found at T'~^^ = 13.99 McV (lower than 
in LGS) and fig^ = 927 A MeV. It is very interesting that 
this reduction of the critical temperatur e ag rees very well 
with the experimental results of Refs. [331 . The further 
properties of the critical point arc given in Table Hvl The 
temperature endpoint is located at T"""^^ = 14.03 MeV 
and A*™^ = 926.6 MeV. We remark that for LGAS the 
temperature endpoint is located on the saturation curve 
(red thick line). This is identical to the location of this 
end-pont for the gas-liquid NCPT in uranium-oxygen 
plasma [6ll - l66j , which is taken as the prototype of NCPT 
for the present study of LGAS (compare Figs. [5] and [5] 
with Fig. 1 of Ref. [Hi). 

In LGAS_fc the two phases are constrained locally to 



0.3. The 



have the same charge fraction Yq 

results arc depicted by the thin gray and black lines in 
Fig. m The two lines also enclose a phase coexistence re- 
gion, which illustrates the non-equality of /is of the two 
phases in the phase coexistence region, due to the locally 
constrained charge fraction (see also the discussion in 
Sec. IV CI and Appendix IC ip . The Gibbs free energy per 
baryon fl is the only chemical potential which is equal in 
the two phases. Furthermore, for isothermal phase tran- 
sitions it is a constant, because the properties of the two 
phases do not change. In contrast to the non-congruent 
phase transition LGAS, for the forced-congruent phase 
transition LG AS_fc /.t^ is dependent on the baryon num- 
ber density ps and given by: 



Pb 



PB 



p%P^a 

PB 



(35) 



which is derived in Appendix IC II 

The phase diagrams as a function of the Gibbs free en- 
ergy per baryon p arc shown in Fig. [5] for both cases 
LGAS and LGAS_fc. The banana-shaped region of 
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FIG. 6. Phase diagram in the plane of temperature and pres- 
sure for cases LGAS and LGAS_fc. In addition to the critical 
point (black dot), the temperature (open square) and pressure 
endpoints (open diamond) are shown in the inlay. Otherwise 
the same notation and depiction as in Fig. |4] is used. 



FIG. 7. The binodal line enclosing the coexistence region in 
the T - pB plane for LGAS and LGAS_fc. Otherwise the same 
notation and depiction as in Fig. |4] is used. 



LGAS is typical for non-congruent liquid-gas like phase 
transitions, see Refs. f6ll - [66j . In Fig. [5] the congruence 
of LGAS_fc becomes obvious. Comparing LGAS and 
LGAS_fc, the phase coexistence region turns into a phase 
coexistence line, when enforcing the local constraint 
for the charge fraction. Furthermore, for LGAS_fc the 
pseudo-critical point (properties listed in Table ITV|l coin- 
cides with the temperature and chemical potential end- 
points. Note that the pseudo-critical point of a forced- 
congruent phase transition obeys Eq. (p4)l . We remark 
that the phase transition line of the forced-congruent 
variant must lie strictly inside the two-phase region of 
the non-congruent phase transition |6ll - [66| . which also 
can be seen as a consequence of Le Chatelier's principle. 
As an exception, both objects could touch each other 
in azeotropic points of the parameter space, as seen for 
LGS. 

Note that for the LGS case, jiB ~ M, since fiq = 0. 
Thus the phase-coexistence line of LGAS_fc in Fig. [5] can 
be directly compared with the one of LGS in Fig. [1] and it 
is found that their shape is very similar. However, states 
on the phase coexistence line of LGAS_fc in Fig. [S] belong 
to two different values of the baryon chemical potential, 
shown by the two thin gray and black lines in Fig. |4l 
If fl in LGAS_fc were changed in a continuous way and 
the phase transition line in Fig. [5] is crossed, jumps 
from the value of the gas phase to the value of the liquid 
phase in Fig. |4l This can be seen as a sign of the enforced 
congruence, in contrast to the azeotropic congruent phase 
transition LGS. 

Figure [S] shows the same scenarios as the previous fig- 
ures, but gives the phase diagrams in the temperature- 



pressure plane. Here it is clear, that the phase transition 
in the forced congruent variant LGAS_fc occurs only at 
a single value of the pressure, which is the same behavior 
as in Fig.[2]for LGS. In contrast, for a given temperature 
in LGAS, there is an extended coexistence range in pres- 
sure, which is enclosed by the SC and BG. The names 
"SG" and "BC" are widely accepted for non-congruent 
evaporation in chemically reacting plasmas (6TI - [65j where 
there meaning is obvious. They arc also most intuitive for 
LGAS for this kind of phase diagram: for a fixed pressure 
of e.g. 10~^ MeV/fm'^ and starting from T = 0, by heat- 
ing the system one will reach the boiling curve, where 
bubbles of gas appear inside the liquid. Conversely, if 
one starts at high temperatures and cools the system 
isobarically, droplets of liquid will form within the gas 
when the saturation curve is reached. In this figure we 
can also identify the pressure endpoint which is located 
on the BC of LGAS. This is identical to the location of 
the pressure endpoint on the BG for gas-liquid NGPT in 
uranium-oxygen plasma [6ll - l66| . which is the prototype 
for our present study of NGPT in LGAS (compare Fig. [6] 
with Fig. 1 in Ref. [Sal)- For LGAS_fc, all three endpoints 
coincide with the critical point defined by Eq. p4p (see 
Fig. 1 in Ref. [H). 



In spite of the similarity of the NGPT in LGAS with 
its terrestrial prototypes [6ll - l66j . the significant differ- 
ence in the topology of P-T diagrams should be stressed 
(compare Fig. [6] above with Fig. 1 in Ref. [Hj). While 
the pressure on both boundaries of the non-congruent 
PT in the uranium-oxygen system (6ll - [66j — boiling and 
saturation curves — tends to zero for the limit T 0, the 
pressure on the boiling curve in NGPT in asymmetric 
[Yq = 0.3) LGPT does not tend to zero when T ^ in 
our case (see Fig.[S]). The same feature was noted already 
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FIG. 8. The charge fraction of the Uquid (dashed hnes) and 
the gas (dotted hnes) for case LGAS as a function of the 
density ps along the binodal hne of Fig. [7] The red hnes 
show the two coexisting phases on the saturation curve of the 
previous figures, the blue lines on the boiling curve, (color 
version online) 



for the NCPT of asymmetric nuclear matter calculated 
with a different EOS (Fig. 3 in Ref. ilM])- The reason for 
this feature is the difference in the physical nature of the 
involved forces which are relevant for the non-con grue nce 
in a chemically-reacting uranium-oxygen plasma |6ll466l | 
and in asymmetric nuclear matter, and also the use of 
Fcrmi-Dirac statistics for the latter. 

In Fig. [7] we show the binodal or phase coexistence 
regions for LGAS and LGAS_fc in the temperature - 
density plane, similar as in Fig. [3] for LGS. Again, the 
non-congruent behavior of LGAS can be identified by 
the non-equivalence of the temperature endpoint and 
the critical point. Conversely, for LGAS_fc the end- 
points and the (pseudo-) critical point coincide. Further- 
more, for isothermal processes of LGAS at temperatures 
rpCP < y < j^TEP so-called retrograde condensation oc- 
curs (see also Ref. QUE]): imagine, e.g., an isothermal 
compression at T = 14 MeV. First one hits the satura- 
tion curve from the left, and a liquid with a larger Yq 
and a larger pB appears inside the gas phase. With in- 
creasing density the volume fraction of the liquid will first 
increase. But for retrograde condensation, for densities 
larger than a certain density, the volume fraction will de- 
crease again, until it returns to zero at the right side of 
the saturation curve. The liquid has disappeared again 
after the phase coexistence region has been crossed. 

The non-congruent behavior of LGAS is further ana- 
lyzed in the following plots of local "chemical" composi- 
tion and density. Fig. [8] shows the charge fractions of the 
two phases which are in coexistence, if one moves along 
the liquid and vapor binodal lines of Fig. [T] In all pre- 
vious plots, the depicted quantities correspond to total 
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FIG. 9. As Fig. [8l but now in dependence of the coexistence 
temperature of the binodal line of Fig. [7| The inlay shows a 
zoom-in to the region around the critical point, including the 
region of so-called retrograde condensation at T > T*-"^. 



thermodynamic quantities, i.e. of the system as a whole. 
In contrast, in Figs. [51 [9l and 1101 we are showing individ- 
ual properties of the two coexisting phases. Now and in 
the following we are using dashed and dotted lines in such 
plots to illustrate this difference. The color coding helps 
to identify the same states in the different diagrams. For 
example in Fig.[8l the ends of the blue curves, which cor- 
respond to T = 0, are given by the highest density of 
coexistence of LGAS in Fig. [7] which is on the blue solid 
line. 

In Fig. [HI one of the two phases always must have Yq = 
0.3, whereas the charge fraction of the second phase is 
not constrained, because its volume fraction is still zero 
on the binodal line. For states on the saturation curve 
one is still in the gas phase, i.e. a = 0, thus Yq = 0.3. 
For states on the boiling curve one is still in the liquid 
phase, a = 1, and Yq = 0.3. The charge fraction is an 
order parameter for LGAS and thus it can be used to 
characterize the two phases, with the identification that 
the gas phase always has a lower charge fraction than the 
liquid, i.e. Yq < Yq . Only at the critical point is equality 
established, Yq ~ Yq. It is interesting to note, that the 
charge fraction of the phase with Yq ^ 0.3 shows the 
same dependence on density before and after the critical 
point. 

In Fig. [S] we also show the charge fractions of the two 
phases along the binodal line, but now as a function of the 
coexistence temperature. By comparing with Figs. [4]- [71 
it is obvious, that for each coexistence temperature there 
are always two points on the binodal line, corresponding 
to two different halves of the binodal line which are sepa- 
rated from each other by the temperature endpoint. For 
each half, two phases with different values of Yq are in 
coexistence. Consequently, in Fig. [SI for each tempera- 
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FIG. 10. The baryon number densities of the liquid (dashed 
lines) and the gas (dotted lines) for case LGAS as a function 
of the coexistence temperature on the binodal line of Fig. [T] 
Otherwise the same notation and depiction as in Fig. [9] is 
used. 
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FIG. 11. Normalized order parameter for chiral symmetry 
restoration/breaking a and order parameter for deconfine- 
ment /confinement "1> versus temperature at zero baryon chem- 
ical potential for symmetric matter of case HIS. 



ture there are always four values of Yq. For T < T , 
one has Yq = 0.3 and Yq > 0.3 on the saturation curve, 

and Y^ < 0.3 and Y^ = 0.3 on the boiling curve. Note 
that the two lines with Yq = 0.3 are on top of each other. 
For <T < T'^^^, both halves belong to the satu- 
ration curve, and thus Yq = 0.3 for both halves, each 
being in coexistence with a different configuration of the 
liquid with > 0.3. 

The previous two figures can be used to identify the 
high degree of isospin distillation of LGAS in the limit 
T — > 0. Let us consider a decompression at T of 
the asymmetric system with Yq = 0.3. Once the boil- 
ing curve is reached vapor bubbles appear which in this 
case consist of pure neutron gas, Yq — 0. Obviously, this 
leads to the distillation of a symmetric liquid by evapora- 
tion of pure neutron bubbles from a boiling asymmetric 
liquid phase. On the other hand, for saturation condi- 
tions ("dewpoint") at T ^ 0, liquid microdrops tend to 
the exactly symmetric composition, Yq — > 0.5. These 
features of the NCPT of case LGAS differ significantly 
from the behavior of the chemical composition (0/U- 
ratio) in NCPTs of uranium-oxygen systems (compare 
Fig. [S] with Fig. 2 of Ref. [U and Fig. 3 of Ref. 

In a similar way as in Fig. [9l in Fig. [10] we show the 
baryon number densities of the two phases for each of 
the two halves of the binodal line as a function of the 
coexistence temperature. Presented in this way, one sees 
that the density is also an order parameter of LGAS, 
whereas the liquid is the phase with the higher density. 
At the critical point, the liquid and the gas have the 
same density and charge fraction and thus cannot be dis- 
tinguished from each other any longer. The discussion of 
this figure is similar as for Fig. [3] For T < T'~^^ there 
are the ps-curves from the liquid (red dashed line) and 



the gas (red dotted line) on the saturation curve, and an- 
other pair of ps-curves from the liquid (dashed blue line) 
and the gas (dotted blue line) on the boiling curve. For 
rpCP < T"^^^ there are still the four different values 
of pb- However, now all points belong to the saturation 
curve. 



B. Deconfinement phase transition 

For the LGPT we used the baryon number density 
and charge fraction as order parameters. For the de- 
confinement and chiral symmetry phase transitions, typ- 
ically the Polyakov-loop $ and the chiral condensate a 
are used, as already discussed in Sec. |IV] The field a 
characterizes chiral symmetry restoration whereas $ can 
be taken as a measure for deconfinement. At finite tem- 
perature and /is = 0, respectively pb = 00 there is no 
first-order phase transition, but a smooth crossover be- 
tween the hadronic (confined, chiral symmetry broken) 
and the quark phase (deconfined, chiral symmetry partly 
restored) [oHi- This is shown in Fig. [11] for HIS, where 
the ratio cr/cro decreases from one to lower values and $ 
goes from zero to a value close to one in a smooth fashion. 
We remark that we define the cross-over temperature 
as the peak of the change of the chiral condensate and 
<f> with T, yielding a value of T'=° = 171 MeV, in accor- 
dance with lattice QCD results [9^. This behavior of the 



We remind the reader that we include anti-particles, therefore 
if T > 0, we have equivalence between pB = and iib = 0. 
Furthermore, for T > 0, iib < corresponds to net anti-matter 
with pB < 0, which is not relevant here. For T = the situation 
is a bit more complicated, because of the LGPT which extends 
down to Pfl = at a constant finite value ol p,B- 
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FIG. 12. Phase diagram in the temperature - baryon chemi- 
cal potential plane for case HIS, for which the deconfinement 
curve (DC) coincides with the confinement curve (CC). The 
black dot marks the critical point. The thick gray line along 
the y-axis shows part of the example trajectory belonging to 
Fig. [Ill 



order parameters corresponds to the example trajectory 
through the phase diagram of HIS shown in Fig. [T^l 

For high enough baryon number densities, the QHPT 
turns into a first-order phase transition. This can be seen 
in Fig. [T2] where we show the first-order phase transition 
line for case HIS, i.e. for heavy-ion collisions of symmetric 
nuclei. The critical point is located at = 165.5 MeV 
and = 383 MeV, again, in accordance with lattice 
QCD results [9^. Its further properties are listed in Ta- 
ble IIVI The topology of this phase transition is the same 
as in LGS, see Fig.[Tl only the typical scales are different. 
For example, the critical point of HIS is at a roughly ten 
times higher temperature. For the QHPT, we will use the 
terms "deconfinement curve" (DC) instead of "saturation 
curve" and "confinement curve" (CC) instead of "boil- 
ing curve" , which we think is more meaningful. If com- 
ing from low densities and temperatures, first droplets 
of denser deconfined quark matter will appear when the 
DC is reached. Conversely, if coming from high densi- 
ties and temperatures, when the CC is reached, the first 
quarks will start to be confined into bubbles of less dense 
hadronic matter. There is an interesting analogy to the 
"ionization boundary curve" and "recombination bound- 
ary curve" of the hypothetical ionization-driven plasma 
phase transition in dense hydrogen at megabar pressure 
range (see, e.g., |l05l - ll07| ). This first-order phase tran- 
sition in neutral hydrogen (H-I-H2) is driven by a jump- 
like ionization (deconfinment) into a predominantly p+e 
plasma state. 

A further general comment about the hadronic and the 
quark phase is required, regarding the assumed full mis- 
cibility of hadrons and quarks, which, e.g., is in contrast 
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FIG. 13. Phase diagram in the temperature - pressure plane 
for case HIS. Otherwise notation and depiction as in Fig. 1121 



to the underlying picture of simple quark-bag models. 
Note that the whole baryon octet and quarks are always 
considered as possible degrees of freedom in the entire 
phase diagram. At sufficiently high temperature, this 
will also lead to the appearance of quarks solutcd in the 
"hadronic sea", i.e., inside what we call the hadronic, re- 
spectively confined phase. On the other hand it is also 
possible that some hadrons survive being soluted in the 
"quark sea", i.e., in the quark or deconfined phase. We 
assume that this inter-penetration of quarks and hadrons 
in the two phases is physical, and it is required to obtain 
the cross-over transition at low baryon chemical poten- 
tial. The two phases are characterized and distinguished 
from each other by their order parameters, which can be 
done in an unambiguous way whenever one has phase 
coexistence. 

HIS is, in principle, a binary system, with baryon num- 
ber and electric charge (respectively isospin) as two glob- 
ally conserved charges (see Table[Tll|). But the phase tran- 
sition in HIS is azcotropic, meaning it is congruent and 
the Maxwell construction can be used, just like symmet- 
ric nucleonic matter in LGS. It is not so obvious as for 
LGS that matter in HIS is an azeotrope, because a whole 
set of particles, including strange ones, is considered (see 
Table IIII[) . However, strange quarks and hyperons do 
not invalidate the relation between isospin symmetry and 
azcotropic behavior, if strangeness is set locally to zero, 
as it is done here. In this case, the total density of strange 
quarks, i.e. in form of unbound strange quarks or bound 
in hyperons, is equal to the total density of anti-strange 
quarks. A more detailed explanation of the matter is 
given in Appendix IB] 

In Fig. [13] the phase diagram in the pressure- 
temperature plane is depicted. Comparing with Fig. [5| 
one realizes an important difference between the QHPT 
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FIG. 14. Phase diagram in the temperature - baryon number 
density plane for case HIS. The deconfinement curve (DC) 
is to the left of the confinement curve (CC). The gray area 
shows the two-phase coexistience region. Otherwise notation 
and depiction as in Fig. 1121 



and the LGPT: the slope of the phase transition line is 
negative. Therefore, the QHPT is not of liquid-gas type. 
With the Clapeyron equation ([5^ one finds that this is 
due to the fact that the hadronic phase, which is less 
dense, has a lower entropy per baryon than the quark 
phase, < s'^ (where we have replaced "/" by "H" 
and "//" by "Q"), which is opposite to the behavior in 
the LGPT. The negative slope of the p — T phase dia- 
gram makes the QHPT fundamentally different from the 
LGPT. This fact (negative slope of the p — T bound- 
ary for QHPT in a symmetric system) is not absolutely 
new (e.g., presentations of I.I. at several conferenceslj, 
L. Satarov, private communication (201 0), J . Randrup, 
private communication (2010), and Ref. |l08j ) but is not 
well-recognized yet. Further investigation and analysis of 
this fundament al di fference between LGPT and QHPT is 
in preparation |l09j . Fig.[T4lshows the coexistence region 
in the temperature - baryon number density plane, in a 
similar way to Fig. [3] for the LGPT. Note that the shape 
of the phase coexistence region of HIS is rather differ- 
ent from LGS. This is once more a manifestation (and 
not the last one) of the fundamental difference between 
LGPT and QHPT. 

In case HIAS of heavy ion collisions with asymmetric 
Goulomb-less matter one has a true binary system. The 
phase transition is non-congruent and the Gibbs' con- 
struction must be used. This is visible in Fig. [TSJ where 
one obtains a phase-coexistence region instead of a phase- 
transition line as in HIS before. We can compare HIAS 
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FIG. 15. Phase diagram in the temperature - baryon chemical 
potential plane for cases HIAS and HIASJc. The slightly 
thicker colored lines belong to HIAS, the slightly thinner gray 
ones to HIAS-fc. The large black dot marks the critical point 
of HIAS and the small black dot the pseudo-critical point of 
HIAS-fc. The filled areas show the coexistence regions. Note 
that for HIASJc the confinement curve (CC) is to the left of 
the deconfinement curve (DC), (color version online) 



with LGAS of Fig. [H Obviously, the phase-coexistence 
region is much narrower than for LGAS, if we compare 
the width in fiB relative to the extension in tempera- 
ture. HIASJc is the forced congruent variant, where the 
charge fraction is constrained locally, Yq = Yq = 0.3, 
so that the Maxwell construction can be used. The light 
gray line is the corresponding DC, and the dark gray line 
the GG, which is partly covering the GG of HIAS. Very 
interestingly, the DG and CC changed order for HIASJc 
compared to HIAS. We will explain this interesting result 
in detail later. 

The phase diagrams as a function of the Gibbs free 
energy per baryon // are shown in Fig. \Wi and as a func- 
tion of pressure in Fig. [ITl Here one sees that HIASJc 
is a congruent phase transition, because the DG and CC 
are identical. Furthermore, the phase transition line of 
HIASJc is inside the phase coexistence region of HIAS, 
as it has to be. Comparing Fig. [16] with Fig. [5] of LGAS, 
again one sees that the shape of the coexistence region 
of HIAS is much narrower. This could be described as a 
weaker non-congruence of the phase transition HIAS. 

Note that in our calculations for HIAS we could not re- 
solve the differences between the temperature and pres- 
sure endpoints and the critical point. In principle, around 
the critical point a similar structure as for LGAS has 
to occur. We predict that the phase coexistence re- 
gions of HIAS around the critical point are smooth and 
two-dimensional, like for LGAS in the inlays of Figs. 0]- 
[71 whereas the temperature endpoint and critical point 
could also have their orders inversed. To be more pre- 
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FIG. 16. Phase diagram in the plane of temperature and 
Gibbs free energy per baryon fi for cases HIAS and HIASJc. 
The thick Ught gray hne gives an example for a trajectory 
through the phase diagram for which we show the order pa- 
rameters in Fig. 1211 Otherwise the same depiction and no- 
tation as in Fig. [15] is used. Note that the phase transition 
line of HIAS-fc is very close to the confinement curve (CC) of 
HIAS. 
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FIG. 17. Phase diagram in the temperature - pressure plane 
for HIAS and HIASJc. The same depiction as in Fig. [15] is 
used. 



else, we expect a round, "banana-shaped" region (see, 
e.g.. Fig. 1 in [HH). However, due to the fact that for 
lower chemical potentials or baryon densities the first- 
order phase transitions HIS and HIAS become rather 
weak, i.e., the phases on both sides of the transition 
become extremely similar, this structure is not as eas- 
ily observed anymore as it is in the LGPT of nuclear 
matter (shown in the previous subsection) and in chemi- 
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FIG. 18. Phase diagram in the plane of temperature and 
baryon number density for cases HIAS and HIASJc. Other- 
wise the same depiction and notation as in Fig. [TS] is used. 
The deconfinement curve (DC) of HIASJc is now to the left 
of the confinement curve (CC). 



cally reacting plasmas [6lH66[ . The highest temperature, 
for which we still could solve the equilibrium conditions 
was at r = 165.3 MeV for = 381 MeV, respectively 
pL = 371 MeV where we used a resolution of roughly 
0.1 MeV in temperature. Because this resolution cannot 
be resolved on the scale of Fig. [121 we take these values 
of T, ^B, and jl as the approximate temperature and 
chemical potentials of the critical point, which are also 
listed in Table IIVI The very narrow phase coexistence 
region shows that the two phases are in an extremely 
similar thermodynamic state. For even higher temper- 
atures, the differences between the two phases become 
smaller than our numerical accuracy, preventing us to 
find a more precise solution. 

This aspect is also very pronounced in Fig. [15] where 
the phase diagrams of HIAS and HIASJc arc shown as 
a function of baryon number density. At T = 0, the DC 
and the CC have very different densities. Conversely, for 
increasing temperatures the extension of the coexistence 
region in density becomes extremely narrow. Note that 
the coexistence region of HIASJc is within HIAS, as it 
must be. 

There is one further important aspect. If we compare 
the phase diagram of the asymmetric system, e.g. Fig. 115) 
with the symmetric system in Fig. 1121 and the corre- 
sponding numbers in Table [TV] one finds that the critical 
points are practically unaffected by the change of the 
asymmetry. T and fiB of HIS and HIAS differ by less 
than 0.5 %. In LGPT, a stronger dependence on the 
asymmetry is observed. For LGS and LGAS the change 
from Yq = 0.5 to Yq = 0.3 leads to a significant shift of 
the critical point, e.g. on the order of 6 % in temperature. 
This is naturally explained by the high temperatures of 
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order 160 MeV around the critical point of HI. The ef- 
fect of the asymmetry becomes extremely weak, because 
the EOS in this regime is dominated by thermal con- 
tributions. Obviously, for fiB — > the EOS would be 
completely independent of the asymmetry. We conclude 
that even though HIAS is in principle a non-congruent 
PT, in our calculations with the Chiral SU(3) EOS model 
this non-congruence is almost negligible close to the crit- 
ical point. This is also visible in Figs. [TSl and [T6l where 
the width of the coexistence region is smaller than the 
line thickness of the curves. 

Let us now come back to the explanation of the in- 
verted ordering of DC and CC in HIAS_fc in Fig. [151 For 
an isothermal compression, corresponding to a horizon- 
tal line through Fig. I18i jl and will increase mono- 
tonically, until the deconfinement curve is reached at a 
certain density ps ~ Pb- This state corresponds to the 
value oi fiB = Mb the hadronic phase on the decon- 
finement curve, i.e. the gray line most to the right in 
Fig. [151 For an isothermal compression, the correspond- 
ing value of fl will remain constant throughout the phase 
transformation, and p = p^ = pf^ . However, inside the 
coexistence region ps will change with density because 
it is given by: 



Pb Pb 



(36) 



with the volume fraction of the quark phase a (see Ap- 
pendix IC ip . For HIAS_fc, the hadronic phase has a 
higher baryon chemical potential than the quark phase, 
p^ > p'^. Therefore an increase of the baryon number 
density inside the phase coexistence region will lead to 
a decrease of pB- For a = 1, i.e. pB ~ Pb, the CC 
in Fig. [T31 is reached. For even higher densities p b will 
increase again. 

Figures [THl and [5D1 show the charge fraction of each 
phase (hadronic and quark) as a function of the temper- 
ature and the baryon density along the phase boundaries. 
For clarity, we again distinguish states on the deconfine- 
ment curve by red color and states on the confinement 
curve by blue. In Figures [TH and [20l one sees that the 
charge fraction in the quark phase is always less than 
or equal to the charge fraction in the hadronic phase, 
Yq < Yq . This is opposite to the liquid-gas PT, where 
the denser phase is more symmetric. 

Note also that Yq can in principle take negative val- 
ues for both the hadronic and the quark phase, due 
to negatively charged hyperons, respectively down and 
strange quarks. In the LGPT we only considered neu- 
trons and protons, so that < Yq < 1. For HIAS we have 
-1 < >Q < 1- Indeed we observe in Fig. [TH that Y^ < 
on the DC for coexistence temperatures below approxi- 
mately 100 MeV. This means that the first quark mat- 
ter droplets which would appear at the deconfinement 
phase boundary of the hadronic phase (DC) are nega- 
tively charged for such temperatures. Remember that on 
the DC only Yq is constrained to the value 0.3, and Yq 
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FIG. 19. The charge fractions of the hadronic [Yq , dotted 
lines) and the quark phase {Yq , dashed lines) as a function 
of the coexistence temperature along the phase boundary, for 
case HIAS. The red lines show states on the deconfinement 
curve DC, and the blue on the confinement curve CC. (color 
version online) 
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FIG. 20. As Fig. 1191 but as a function of the coexistence 
baryon density. 



is set by the equilibrium conditions. Conversely, if the 
phase transition is crossed coming from the high density 
side we have Yq ~ 0.3 on the CC. The hadronic phase on 
the CC shows the opposite tendency than quark matter. 
For low temperatures it approaches a rather symmetric 
configuration with Yq close to 0. 5. Similar features have 
been discussed, e.g., in Ref. |llO| for a simple quark-bag 
model. Fig. [20l appears quite complex and looks different 
to the equivalent Fig. [51 of LGAS. This can be explained 
by the non-monotonous behavior of the density as a func- 
tion of the coexistence temperature in case HIAS, visible 
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FIG. 21. Normalized order parameter for chiral symmetry 
restoration/breaking a (black lines) and order parameter for 
deconfinement/confinement $ (red lines) versus the Gibbs 
free energy per baryon for case HIAS. Inside the phase co- 
existence region, the values of the order parameters in the 
quark phase (dashed lines) are different from the ones in the 
hadronic phase (dotted), and change during the phase trans- 
formation. HIASJc is the forced congruent variant, shown 
by solid lines. All order parameters are calculated for the 
trajectory through the phase diagrams shown in Fig. [16] with 
T = 0.15/i. (color version online) 



in Fig. [THl 

The two figures [12] and [201 also show how the charge 
fraction in both phases goes to 0.3 when approaching 
the critical point. At the critical point, the two phases 
are identical, i.e. have the same charge fraction, density, 
scalar field $, chiral condensate, etc. To better under- 
stand the dynamics of such a non-congruent phase tran- 
sition, in Fig. [16] we included an example for a trajec- 
tory through the phase diagram for which we show the 
order parameters in Fig. [^ The trajectory was cho- 
sen to follow T ~ 0.15/2 so that the phase coexistence 
region is crossed at temperatures around 900 MeV and 
/t ~ 130 MeV. This trajectory could e.g. be realized dur- 
ing the decompression of the quark-gluon plasma in a 
heavy-ion collision. We show the behavior of the order 
parameters for HIAS and HIAS_fc to illustrate the dif- 
ferences between a non-congruent and congruent phase 
transition. For HIASJc, when the phase transition line 
is crossed, there are jumps in a and as expected for a 
first-order phase transition (compare also with Fig. [TT]) . 
The quark phase is characterized by having a larger value 
of $, showing that it represents the deconfined state. As 
can be seen from Eqs. ([1]) and ([2]), this increases the ef- 
fective mass of baryons, but decreases the effective mass 
of quarks. Note that a, which is increasing with decreas- 
ing /i in both of the two phases, locally decreases at the 
phase transition (going from right to the left). This ef- 
fect comes from the fact that both order parameters are 
connected through the effective masses of the particles. 
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FIG. 22. Phase diagram in the temperature - baryon chemical 
potential plane for cases NSLCN and NSGCN. The horizontal 
thick gray line at T = is a trajectory through the phase 
diagram, for which we show the order parameters in Fig. 1231 
The large black dot is the critical point of NSGCN, the small 
one the pseudo-critical point of NSLCN. (color version online) 



Also the baryonic density of the quarks (not counting the 
contribution from U) is less than the one of the hadrons. 
This leads to a decrease of the scalar field across the 
phase transition. 

In Fig. [21] it can be seen that in the non-congruent 
phase transition HIAS the behavior of the order param- 
eters is more complex. Within the phase coexistence re- 
gion, i.e. within the DC and CC shown in Fig. [TB] the 
hadronic and the quark phases are in coexistence. The 
two phases are spatially separated and in each of the 
two phases one has different values of the order param- 
eters. At the onset of the phase transition at the con- 
finement curve one has a = 1, i.e. the volume fraction 
of the quark phase is still one and the volume fraction 
of the hadronic phase is 0. Inside the phase coexistence 
region, the quark phase has a larger value of but also 
a slightly increased chiral condensate a, as observed for 
HIAS_fc before. With decreasing /i, not only the vol- 
ume fraction of the quark phase decreases, but also the 
properties of the two phases change. One clearly sees 
that chiral symmetry breaking proceeds in both of the 
two phases with decreasing density. On the other hand 
is decreasing, and 3>*^ is slightly increasing. When 
the deconfinement curve is reached, a = and only the 
hadronic phase is left. After this in Fig. [211 one sees that 
<I>-^ and become equal to values obtained for HIAS_fc. 

Finally, we discuss results for the quark-hadron phase 
transition in neutron stars shown in Fig. [5^ Let us re- 
iterate that in this paper we consider the quark-hadron 
PT in multi-component matter of neutron stars of macro- 
scopic coexisting phases (i) with local charge neutrality 
and (ii) in the variant of global charge neutrality within 
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the Coulomb-less approximation. It is just the Coulomb- 
less approximation which justifies the use of Gibbs con- 
ditions instead of Gibbs-Guggenheim conditions, which 
are valid for true macroscopic Coulomb systems (e.g., 
Ref. [1^). Local and global charge neutrality constraints 
lead to forced-congruent ( "Maxwell" ) and non-congruent 
("Gibbs") phase transitions. As discussed in Sec. Ill Al 
the choice between local or global charge neutrality can 
be associated with the unknown surface tension between 
the two phases. 

The phase diagram of neutron stars looks similar to 
the one shown in Fig. 1161 In contrast, the dimensionality 
of HIASJc in Fig. [15] is different to the one of NSLCN in 
Fig. [211 This is easy to understand, because = fig 
is vahd for NSLCN and NSGCN, and also for HIAS, but 
not for HIAS_fc. The critical points of the neutron star 
cases are at lower /is and slightly higher T than in cases 
HIAS and HIS, i.e., the phase transition regions extend 
to higher temperatures. Because asymmetry and charge 
neutrality do not seem to have a major effect at such 
high temperatures, this is related to the different treat- 
ment of strangeness (see Table IIIip . At T = the width 
of the phase coexistence region of HIAS is 52 MeV, for 
NSGCN it is 70.2 MeV. With the chiral SU(3) model 
used here, the hadron-quark phase transition and the 
non-congruent features at T ~ seem to be stronger for 
neutron star matter than for matter in heavy-ion colli- 
sions. This is a result of the larger asymmetries obtained 
in beta-equilibrium due to the large degeneracy of elec- 
trons, and also because of the different constraints used 
for strangeness. 

The critical point of NSGCN is at approximately 
T^P = 168.9 MeV and = 224 McV, the one of 

NSLCN at T'^P = 168.6 MeV and ^g^ 247 MeV (see 
also Table lIVp . Interestingly, they are only slightly dif- 
ferent. It shows that the treatment of electric charge 
neutrality plays only a minor role for the location of crit- 
ical points of the QHPT at the typical high temperatures. 
The local constraint applied in HIAS_fc which does not 
allow isospin diffusion, has a slightly larger effect at high 
temperatures. 

The phase transition line for NSLCN with enforced 
local charge neutrality must lie strictly within the phase 
coexistence region of NSGCN, in accordance with general 
rules for non-congruent phase transitions (e.g., Refs. [611 - 
[66j). Both boundaries can touch each other in points 
of azeotropic composition. Here this would mean that 
the two phases were charge neutral a priori. Actually 
this is not the case, but they only touch because of the 
thickness of the lines in the figure. For NSLCN, fiq will 
not behave continuously when the phase transition line 
is crossed, as already discussed at the end of Sec. IV C| 
because /Xg is different from fig. This is in accordance 
with general properties of phase coexistence of charge 
neutral, macrosco pic p hases in C oulomb systems (see, 
e.g., Refs. [5^ llll| and |ll2| ): any phase-interface 
in macroscopic equilibrium Coulomb systems is accom- 
panied by a finite difference in the average electrostatic 
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FIG. 23. As Fig. [m but now for cases NSLCN and NSGCN 
of the example trajectory at T = 0. 



potentials of both coexisting phases (Galvani pot ential), 
see also Sec. Ill Al and compare also with Ref. jll3| . How- 
ever, the total charge chemical potential /ig still can be 
related to the local chemical potentials fiQ and /Iq , if it 
is seen as a function of baryon number density. Based 
on two different assumptions about the implementation 
of local charge neutrality, in Appendix IC 31 we derive the 
following expressions: 



MQ = 



Pb 



tor PQ = Pq 
for Y" = 



(37) 



It would be more intuitive to assume the first of the two 
conditions. These expressions can be used to determine 
the total chemical potential of charged particles inside 
the two-phase mixture. 

Fig. [^ also includes an example trajectory through the 
phase transition region at T = 0. This trajectory could 
correspond e.g. to the spatial structure inside a neutron 
star, or better called hybrid star. In Fig. [23] we show 
the behavior of the order parameters for cases NSLCN 
and NSGCN along the trajectory. For T = 0, we find 
<I>'^ = 0, i.e. the hadronic phase is completely confined. 
For finite temperature this is different as can be seen in 
Figs. [TT] and [m where one has some sequential deconfine- 
ment already in the hadronic phase. Otherwise similar 
features are observed as in Fig. [2TJ The obtained jumps 
show that we have a first-order phase transition. The dif- 
ferent values of the order parameters distinguish the two 
spatially separated phases. If the phase transition is non- 
congruent, one obtains a phase coexistence region, and 
the phases change their properties, including the order 
parameters, during the phase transformation. 
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TABLE IV. Approximate location of the critical and pseudo- 
critical points of the different scenarios. 

case r [MeV] pg [fm'^'] fig [MeV] p [MeV/fm^] 



LGS 


14.75 


0.046 


912.4 


0.205 


LGAS 


13.99 


0.049 


927.4 


0.241 


LGASJc 


12.68 


0.044 


928.3 


0.171 


HIS 


165.5 


0.24 


383 


73.1 


HIAS 


165.3 


0.22 


381 


68.6 


HIASJc 


156.5 


0.44 


611 


95.8 


NSLCN 


168.6 


0.16 


247 


162.8 


NSGCN 


168.9 


0.14 


224 


161.1 



VII. SUMMARY AND CONCLUSIONS 

In this article we investigated the hquid-gas phase 
transition (LGPT) described by the FSUgold relativistic 
mean-field model and the hadron-quark (HQ) or decon- 
finement phase transition with the Chiral SU(3) model. 
We did not take into account any explicit Coulomb in- 
teractions or other finite-size effects, but always consid- 
ered phase coexistence of macroscopic phases. Different 
physical systems were investigated: heavy ion collisions 
of symmetric nuclei at low (LGS) and high collision en- 
ergies (HIS) and of asymmetric nuclei with Z/A = 0.3 
at low (LGAS) and high collision energies (HIAS). Fur- 
thermore, we also studied the QHPT in neutron stars 
(NS). The main goal of our work is to characterize the 
first-order phase transitions in the different systems re- 
garding their possible non-congruence and to identify the 
typical non-congruent aspects. This characterization and 
classification, and our qualitative results should not de- 
pend on the model used for the calculations, but should 
be valid in a rather general way. In our study of non- 
congruent phase transitions we used essential features of 
such type of phase transitions obtained from terrestrial 
applications with high-temperature chemically reacting 
plasmas [6TI - [66j . 

A non-congruent, first-order phase transition occurs 
for a non-unary system with more than one (globally) 
conserved charge where the charge ratios change dur- 
ing a phase transformation. In principle, all the inves- 
tigated scenarios contain several conserved charges, with 
at least baryon number and electric charge. However, 
due to isospin symmetry, for the symmetric case (LGS 
and HIS), matter always stays symmetric, even if there is 
phase-coexistence, which is called "azeotropic" behavior. 
Only these two cases lead to congruent phase transitions. 
To further illustrate the differences between a congruent 
and a non-congruent phase transition, we also consid- 
ered the cases LGAS_fc and HIASJc, where we forced 
the system to have a congruent phase transition, by con- 
straining the local charge fractions to be equal in the 
two phases. From a physical point of view the forced- 
congruent regime of phase transformations is not com- 



TABLE V. Classification of the equilibrium conditions and 
type of first-order phase transition in the different scenarios. 



case 


equilibrium conditions 


type of PT 


LGS 


Gibbs/Maxwell (equiv.) 


congruent (azeotrope) 


LGAS 


Gibbs 


non-congruent 


LGASJc 


Maxwell 


forced congruent 


HIS 


Gibbs/Maxwell (equiv.) 


congruent (azeotrope) 


HIAS 


Gibbs 


non-congruent 


HIASJc 


Maxwell 


forced congruent 


NSLCN 


Maxwell 


forced congruent 


NSGCN 


Gibbs 


non-congruent 



pletely artificial, but corresponds to the so-called "frozen 
diffusion approximation" . In a similar way we consid- 
ered enforced local charge neutrality in neutron stars in 
case NSLCN instead of global charge neutrality NSGCN. 
In Table |V] we summarize the classification of the phase 
transitions in the different physical systems. 

In the astrophysics community, congruent and non- 
congruent phase transitions are usually called "Maxwell" 
and "Gibbs" phase transitions, referring to the way in 
which phase coexistence is calculated. Because of the 
constant charge ratios and the fact that all proper- 
ties of the two phases do not change for an isother- 
mal phase transformation, a simple Maxwell construc- 
tion can be used for the congruent case. Conversely, for 
non-congruent phase transitions the more complicated 
Gibbs construction is necessary to solve the set of ther- 
mal, mechanical and multiple chemical equilibrium con- 
ditions. Due to their prominent role in neutron star 
physics, NSLCN and NSGCN are considered as the two 
most representative scenarios. However, as we demon- 
strated, they are only examples illustrating the further 
classification of first-order phase transitions to be either 
congruent or non-congruent. 

The distinction between congruent and non-congruent 
phase transitions is crucial for the dimensionality of phase 
diagrams. For a given temperature T, a congruent phase 
transition occurs at a single value of the Gibbs free energy 
per baryon /i. Therefore in the T-/i-plane, one obtains 
a phase transition line, as well as in the T-P- and P- 
/i-planes, and generally in any pair of intensive thermo- 
dynamic variables [ssj . Conversely, for a non-congruent 
phase transition calculated with the Gibbs construction 
and a given chemical composition, in any pair of intensive 
thermodynamic variables one obtains two-dimensional 
"banana-like" phase-coexistence regions (see, e.g.. Fig. 1 
in Ref. [IHl)- One of the conclusions from the present 
study is that all asymmetric systems are non-congruent, 
if no additional local constraints are enforced. The corre- 
sponding phase diagrams should always be regions, and 
not lines. 

Inside the phase coexistence region, the order param- 
eters (e.g. density and charge fraction for LG, density. 
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charge fraction, chiral condensate, and $-ficld represent- 
ing the Polyakov-loop for QHPT) have different values in 
the two phases, which is the actual definition of a first- 
order phase transition. We demonstrated the behavior 
of some of these order parameters on the binodal bound- 
aries and within two-phase regions, and how they become 
equal at the critical point. Furthermore, for the QHPT 
this was also illustrated for the chiral condensate and the 
<i>-fLeld for trajectories through the phase diagram in a 
high energy heavy-ion collision and in a neutron star. For 
a congruent phase transition, the different values appear 
as a jump when crossing the phase transition line and 
moving from one phase to the other. For non-congruent 
phase transitions, at the phase coexistence boundary the 
second phase appears, and its volume fraction increases 
until the first phase has disappeared completely when 
the phase coexistence region has been crossed. During 
such a non-congruent phase transformation, both phases 
continuously change their thermodynamic properties and 
order parameters. 

Local charge neutrality for macroscopic phase coex- 
istence in neutron stars under Maxwell conditions has 
been applied in a number of published works. It is well 
known, that the local charge chemical potentials behave 
discontinuously in such a transition, while the general- 
ized non-local electro-chemical potentials are continuous 
in accordance with generalized Gibbs-Guggenheim con- 
ditions (e.g. Ref. 15a])- The discontinuity is associated 
with the so-called Galvani potential of the phase interface 
in Coulomb systems (e.g., Refs. 65 1 and jll4j ). In the 
present article, for the Coulombless approximation we de- 
rived expressions for the total charge chemical potential 
inside the phase coexistence region as a function of 
the total baryon number density, which behaves continu- 
ously. The new expressions show explicitly that the total 
chemical potentials on the binodal surface correspond to 
the ones of the dominating phase if its volume fraction is 
unity. We also derived the total chemical potentials for 
the other phase transitions where local constraints were 
applied. 

Another fundamental aspect of non-congruent phase 
transitions is the location and properties of the criti- 
cal point. For congruent phase transitions, the criti- 
cal point is located at the end of the phase transition 
line. Conversely, the phase coexistence region of a non- 
congruent phase transition has several (topological) end- 
points: namely for temperature, pressure and chemical 
potential, which are all different (e.g., Ref. [1^). The 
critical point, as the point of thermodynamic singularity, 
and where the two phases arc identical, docs not coincide 
with these endpoints. For the LGPT we found that the 
temperature endpoint is located on the saturation curve, 
i.e. the critical point is located at higher densities than 
the temperature endpoint, while the pressure endpoint is 
located on the boiling curve. This is the same arrange- 
ment as for non-congruent phase transitions in terrestrial 
chemically reacting plasma (e.g., Fig. 1 in Ref. [l^). The 
properties of all (pseudo-) critical points found in the 



present study arc listed in Table IIVI 

For the QHPT we could not resolve the endpoints and 
the critical points because we reached the limit of our 
numerical accuracy. But at the same time, this is related 
to another important finding from our study: at low tem- 
peratures, the asymmetry has a significant impact on the 
phase diagrams of the QHPT. Conversely, for high tem- 
peratures, close to the critical point, this impact is van- 
ishingly small because the EOS is dominated by thermal 
contributions. As a consequence, the non-congruent fea- 
tures of HIAS are almost unnoticable close to the critical 
point. This is in contrast to LGAS of nuclear matter 
considered here and LGPT in terrestrial chemically re- 
acting plasma [6ll - l66j for which the non-congruence is 
significant for all conditions. 

Another important difference between LGPT and 
QHPT comes from the phase diagrams in the 
temperature-pressure plane. For ordinary Van-der- 
Waals-like phase transitions of liquid-gas type, and also 
for the numerous variants of LGPT in nucleonic matter 
studied in other papers, the slope of the P-T-phase tran- 
sition line is positive, which we also confirm in our study. 
Conversely, for the QHPT we found that the slope is neg- 
ative, which leads to the conclusion that this phase tran- 
sition is not of liquid-gas type. The Clapeyron-equation 
shows that this peculiar feature is a result of the gener- 
ally higher entropy per baryon of the quark phase than 
of the hadronic phase during phase coexistence. We also 
found that the denser phase of the QHPT (the quark 
phase) is more asymmetric than the more dilute phase 
(the hadronic phase), which is also opposite to the be- 
havior in the LGPT. 

In the future, we could extend our work to study phase 
diagrams of proto-neutron stars and supernovae. Those 
objects are special, not only for having higher temper- 
ature in comparison to neutron stars, but also for hav- 
ing a higher fixed lepton content due to the presence of 
trapped neutrinos. Such features transform an already 
binary phase transition into a ternary one due to the in- 
troduction of the lepton number conservation, with po- 
tentially interesting results. It was already shown that a 
phase transition to deconfined matter prio r to the super- 
nova explosion can have observable effects jll5{ . Another 
aspect which we did not discuss in detail are the con- 
sequences of different constraints regarding strangeness 
conservation and possible strangeness distillation. 

It would be interesting to explore further the experi- 
mental consequences of non-congruence in heavy-ion col- 
lisions. In our investigations we have found that the non- 
congruence of the QHPT increases with decreasing tem- 
perature. Therefore, the possible non-congruent features 
could become particularly relevant for the heavy-ion ex- 
periments at the future FAIR facility at Darmstadt or 
NICA at Dubna, which both aim to probe asymmetric 
matter at high densities and low temperatures. 
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Appendix A: Critical points in a "two-EOS" 
description 

It is well-known (e.g., Ref. [H) that in a "two-EOS" 
description of PTs, where two different EOS models 
are used for the two phases in coexistence, there can- 
not be a termination point or critical-point of a phase 
transition line. This is due to the standard argument 
of the impossibility of a continuous and smooth trans- 
formation from one phase into another in a two-EOS 
approach (e.g., commonly used for crystal-fluid coexis- 
tence). This is in contra st to the claims in the recent 
studies of, e.g., Refs. (33l [TT6l - [TT8 |. The previous state- 
ment can be supported with the following geometrical 
argumentation. For symmetric matter, the two phases 
give two different planes in the parameter-space of pres- 
sure P, temperature T and baryon chemical potential 
/ifi. Where the two planes intersect, one has phase coex- 
istence. This intersection of two two-dimensional planes 
in a three-dimensional parameter space either has to be 
a closed or an infinite curve. A termination point of the 
intersection curve is impossible, as long as the planes do 
not show any discontinuities, which would be rather un- 
physical. This argumentation can also be generalized to 
asymmetric matter, where the charge chemical potential 
/ig appears as another dimension. 



Appendix B: Azeotropic behavior for symmetric 
matter in HIS and LGS 



The azeotropic behavior of HIS and LGS can be ex- 
plained result of isospin symmetry in the follow- 
ing way. First of all, it is important to note that from 
Yq = 0.5 it follows that the total third component of 
the isospin is zero, /a = 0, if there is no net strangeness, 
S" = 0. Furthermore, the considered particles (neutrons 
and protons for LGS, u-, d-, s-quarks and the baryon 



octet for HIS) can be grouped in pairs of isospin partners, 
where the isospin partners have identical baryon number, 
strangeness and mass, but opposite third component of 
the isospin, and some remaining particles with zero third 
component of the isospin. If, in addition, the interactions 
arc built in an isospin symmetric way, namely that they 
arc identical for the isospin partners at /a = 0, it follows 
that /a = if and only if the isospin chemical potential 
is also zero, /i/, — ^ — 0. For the constraints used 
in LGS and HIS all these conditions are fulfilled by the 
two applied EOS models FSUgold and Chiral SU(3). On 
the other hand, it can be shown that the charge chem- 
ical potential is equal to the isospin chemical potential, 
M/3 = I^Q- Thus, we obtain that = 0, if Yq = 0.5 and 
5 = 0, but independent of temperature and density. For 
two-phase coexistence in HIS and LGS we use the con- 
straints = S'^ — and Yq = 0.5. Therefore we get 
/iQ = = and Yq = Yq^ = 0.5, which shows that the 
two phases remain symmetric and that the phase transi- 
tion is therefore azeotropic. 



Appendix C: Total chemical potentials inside the 
phase-coexistence regions 

In all HI and LG cases local constraints are applied 
for two-phase coexistence. Therefore, it is necessary to 
refine the definition of the total chemical potentials inside 
the two-phase mixture, and to relate them to the local 
chemical potentials already defined in Table HH 



1. HIAS_fc and LGAS.fc 

Let us start with cases LGASJc and HIAS_fc. The 
constraints listed in Table IIIII are equivalent to consid- 
ering B = const., Q = const. = YqB, Yq = Yq^ , and 
= S^^ = 0. Note that this set of constraints also 
gives Yq — Yq ^ Yq^. The constraint for zero local 
strangeness, = S^^ = 0, is formulated with local ex- 
tensive variables. However, it would be equivalent to 
consider Yg = Yg^ and S* = 0, or = p^f and S = 0, 
where we define Yg ~ /B^ , pg ~ S-^ /V-' , and anal- 
ogous expressions for phase //. In conclusion, we have 
the simple conservation of the total charges, which are 
fixed to some values, plus two additional local constraints 
Y^ = Y^^ and F/ = Yf, or Y^ = Y^^ and p'g = p^/. 
Let us first use the first of the two formulations of the 
two local constraints. 

Now we can define the total baryon chemical potential 
inside the two-phase mixture as the following derivative 



Mb 



dF 
dB 



T.V,S.Q.Y'=Y'J,Y'=Y, 



(CI) 



This is nothing but the definition of /i^ outside the two- 
phase mixture given in Table |T1 but with the additional 
local constraints taken into account. This is what we 
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mean by the total chemical potentials inside the two- 
phase mixtures. For the definition of /xg, "Q" and "B" 
have to be exchanged in Eq. (jCip . for fxs, "5"' and "i?". 
For p., only "Q" has to be replaced by "Yq" in the list 
of constant variables. Note that this definition of /i gives 
the same relation 



(C2) 



(C3) 



obtained previously for single phases. Let us introduce 
X as an abbrcvation for the set of variables kept constant 
and the local constraints 

OF 

For the definition of the total chemical potentials of other 
cases than HIAS_fc and LGAS_fc, one only has to ex- 
change the local constraints in X. 

Now we want to evaluate the above expression for /is 
in HIAS_fc and LGAS_fc. Inside the two-phase mixture, 
the total free energy F is given as the sum of the two 
phases 



F^\T,V 



II j^ii 



,5" g"), 

(C4) 

where we used thermal equilibrium, Eq. (fT2|l . With the 
chain rule and the definitions of the local chemical po- 
tentials of Table we obtain 



dB 
dB' 



dV 



II 



X 



dB 
'dB 



dB 



dB 
dB" 



II 



dB 

ds" 



p 



dB 
dQ" 



X 



dB 



(C5) 



Because of pressure equilibrium, Eq. the first two 

terms sum up to zero. Next we use = YgB^ for the 

dS' 

X 



expression 



dS' 
~dB 



= Y, 



J dB' 



dB 



B 



idYi 



dB 



= B 



I dYi 



dB 



,(C6) 



where the second equality comes from Yg = 0. If we have 
Ys = Y^', we also have = y/^ = Ys, and thus 



B 



I dYi 



dB 



B 



I &Y^ 
dB 



B^ 
B 



-Ys 



(C7) 



To make use of the local constraint Yq = Yq' in X, 
we replace Q' by YqB' , and Q'' by Yq'B'' , and use 
Yq^Y'^ Y". This gives 



dB' 



dB 
dB" 



X 



dB 



II .,II\ 



X 



- —Y'n' - —Y"u" 



(C9) 



Now we can use the inter-phase equilibrium condi- 
tions (^5]) and (P^ . and the definitions of jl' and jl" , 
(Eq. dnil)), to obtain 



B^ 
" B 
B' 



B^ 
B 



Ms 



B 
~B 

B" 



B'^ 
3 
B^ 
B 



B" 
B 



Y"ii" 



' B 

This can be written as 



MB = (l-a)^/4 

PB 



PB 



(CIO) 



(Cll) 



where a is the volume fraction of phase // in the two- 
phase mixture 



a = 



PB 



p'b 



P'J 



Pb 



(C12) 



For Pb = Pb one gets pB ^ Pb^ ^-^id for pB = p" Pb = 
Pb , i.e. the correct limits are obtained. Furthermore, if 
Pb = Pb '^^^ finds that pB = Pb ~ Pb^ independent of 
Pb- With a similar derivation one obtains the following 
expression for the total charge chemical potential 



(C13) 



It turns out that one obtains the same expressions (jClip 
and (jClSp . if one uses the strangeness constraint in the 
form pg — p'g instead of Yg = Y" . Furthermore, note 
that Eqs. ()Clip . (jC13j) . and (fC2|) also allow one to ex- 
press p in terms of the local chemical potentials. 

Interestingly, for the total strange chemical potential 
the form of the strangeness constraint makes a difference 



In conclusion, we obtain 



also ^ 



X 



Pb 



ds' 

dB 



0, and in the same way 



= 0. Thus we are left with 



dB' 



dB 
dQ' 



dB 



Pb 



X 



Pq 



dB 



II 



dB 
dQ" 



dB 



p'J 



X 



P'^ 



(C8) 



_ J (1 - a)p'g + ap'J, tor py = pg 

'''{{1- «)ff + "f//, for Yi . r/^'"^ 

Note that the two forms give identical results if a = or 
if a = 1, i.e., at the phase boundaries. 

Eqs. (jClip . (jC13|) . and (jC14p can be used to determine 
the chemical potential of particles for the two-phase mix- 
ture as a whole. For example, for protons and neutrons 
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one obtains using Eq. ([7]) 



Appendix D: Gibbs free energy 



Pb Pb ' 

PB PB 



(C15) 
(C16) 



For A's the constraint Yg = would, e.g., lead to 

MA = /^i + a— Ma ■ (C17) 

Pb Pb 



2. LGS, LGAS, HIS, and HIAS 

In cases LGS, LGAS, HIS, and HIAS the only local 
constraint is = S^^ = 0. Here one finds 



Pb — Pb — Pb j 
Pq = Pq ^ Pq i 



(C18) 
(C19) 



PS = 



for pi; = p{;' 



{1 - a)pl; + aplJ , 



(C20) 



3. NSLCN and NSGCN 

For case NSGCN, which uses the Coulomb-less approx- 
imation and leads to a non-eongruent FT, all three local 
chemical potentials have equal values in the two phases 
and one obtains 



Pb - Pb - Pb i 
P-Q ^ Pq = Pq > 



Ms 



^pi = pi'{^ 0). 



(C21) 
(C22) 

(C23) 



In Table IIIII local charge neutrality constraint in 



NSLCN was formulated as Yq 



Yq = 0, i.e. in terms 



of the local electric charge fractions. Here one has a 
similar situation as for the local strangeness constraint 
in Appendix IC II it would be equivalent to consider 
Pg = p^ = instead of Yq = Yq = 0. After a sim- 
ilar derivation as in Appendix IC II one obtains 



Pb 


= Pb = 


PS 


= Ps = 






PQ 


-It- 



p'J 



a)pl 



(C24) 
(C25) 



for Pq 
for 



^Q - 



Pq 

^Q 



(C26) 



For the determination of pQ it makes a difference whether 
one assumes equal charge fractions or equal charge den- 
sities in the two phases. It would be more intuitive to 
assume the first of the two conditions, namely that ad- 
ditional charge is distributed uniformly, because here we 
have in mind that the Coulomb forces are the reason for 
local charge neutrality. 



Let us first assume we are outside of the two-phase 
coexistence region, i.e., there is only one phase. The 
Gibbs free energy G can be obtained via a Legendre- 
transformation from the (Helmholtz) free energy F 



G=F+PV 



(Dl) 



If G = G(T, P, {Ni}) is seen as a function of the particle 
numbers Ni, which are fixed by Eq. ([7]), one obtains 



G = ^ N,p, 

i 

= BpB + Qpq + Sps 



(D2) 
(D3) 



This shows that G can also be seen as a function oi B, Q, 
and S, i.e., G = {T, P, B,Q, S). Furthermore, because 
we consider either S — ot ps — 0, we obtain from 
Eqs. dnH, O and © 



G = BpB + Qpq 
= Bp , 



(D4) 
(D5) 



which shows that p is indeed the Gibbs free energy per 
baryon, respectively the two definitions are equivalent. 
Inside the two-phase mixture we have 

G = G(T,P,{A^/},{7V/0) 

= G\T, P, {nI}) + G"{T, P, {iV/^}) 

= 5]iV/A.f +^iV/Vf' (D6) 



B'p'b + 
+S'p's 



B"p'J ^ 



Q'pq + Q 



Pq 



(D7) 



where we used Eqs. ((TH)) . Interestingly, the total chemical 
potentials given in Appendix [C] also fulfill the following 
equalities 



G 



(D8) 
(D9) 



BpB + Qpq + Sps 
^G{T,P,B,Q,S) . 

Note that 

G^(T, P, {Nl }) = bVb + S'p's + Q'pq , (DIO) 

and because we always consider either 5*^ — S^^ = or 
Ps = p¥ = 0, we obtain from Eq. (j3ip 

(Dll) 
(D12) 



G' - sVb + Q'pq 



B'p' 



so the local p^ is also equal to the local Gibbs free energy 
per baryon. 

We remark that the Gibbs free energy per baryon is 
important because it is just this quantity which must be 
equal in the case of forced-congruent phase coexistence 
(under Maxwell conditions), in analogy to the specific 
Gibbs free energy used for terrestrial chemically reacting 
plasma. 
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